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ABSTRACT 

We have carried out an in-depth study of low-mass X-ray binaries (LMXBs) detected in the nearby 
lenticular galaxy NGC 3115, using the Megasecond Chandra X-Ray Visionary Project observation 
(total exposure time 1.1 Ms). In total we found 136 candidate LMXBs in the field and 49 in globular 
clusters (GCs) above 2 a detection, with 0.3-8 keV luminosity Lx ~10 36 -10 39 erg s _1 . Other than 
13 transient candidates, the sources overall have less long-term variability at higher luminosity, at 
least at Lx > 2 x 10 37 erg s _1 . In order to identify the nature and spectral state of our sources, we 
compared their collective spectral properties based on single-component models (a simple power law 
or a multicolor disk) with the spectral evolution seen in representative Galactic LMXBs. We found 
that in the Lx versus photon index Tpl and Lx versus disk temperature /cTmcd plots, most of our 
sources fall on a narrow track in which the spectral shape hardens with increasing luminosity below 
Lx ~ 7 x 10 37 erg s _1 but is relatively constant (Ppl~ 1-5 or /cTmcd~L 5 keV) above this luminosity, 
similar to the spectral evolution of Galactic neutron star (NS) LMXBs in the soft state in the Chandra 
bandpass. Therefore we identified the track as the NS LMXB soft-state track and suggested sources 
with Lx < 7 x 10 37 erg s~ 4 as atolls in the soft state and those with Lx > 7 x 10 37 erg s _1 as Z 
sources. Ten other sources (five are transients) displayed significantly softer spectra and are probably 
black hole X-ray binaries in the thermal state. One of them (persistent) is in a metal-poor GC. 

Subject headinqs: X-rays: binaries — globular clusters: general — Galaxy:stellar content — X-rays: 
individual (NGC 3115) 


1. INTRODUCTION 

In a low-mass X-ray binary (LMXB), a neutron star 
(NS) or a stellar-mass black hole (BH) accretes matter 
from a Roclie-lobe filling, low-mass companion star via 
an accretion disk. Co p ious X-rays a re produced in the 
inner disk (lSliakura fe Sunvaevlll973l) . In the case of NS 
LMXBs, there is also strong X-ray emission from the 
boundary layer formed by the set tling of the accretio n 
flow onto the NS s u rface dlnogamov fc Sunvaevl 119991 
iPopham fe Sunvaevl 1200 lH . Our knowledge of LMXBs 
was revolutionized thanks to 16 years of intensive obser¬ 
vations of such objects in our Gal axy by the Rossi X-ray 
Timing Explorer (RXTE. iBradt et al.lll993l l. 

BH LMXBs constitute the majority of BH X-ray bi¬ 
naries (BHBs) known in our Galaxy, though a few 
BHBs with a high-mass stellar companion also ex¬ 
ist. BHBs exhibit three main X-ray spectral states: 
the hard state, the thermal state, and the steep 
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power law (SPL1 state (iRemillard fc McClintockl 120061 

iMcClintock &; Remillardl12006T ) . These states are nor¬ 
mally described with spectral models consisting of two 
main continuum components: a stan dard thermal multi¬ 
color disk (MCD, diskbb in XSPEC. lArnaudl 11996 1 and 
a Comptonized component (often modeled with a single 
power law, PL, or the Comptonization model comptt in 
XSPEC). The hard state, which tends to occur below 
several percent of the Eddington luminosity (Le<m)j is 
characterized by a strong Comptonized component dom¬ 
inating the spectra at least above ~2 keV and extend¬ 
ing to more than 100 keV. A weak, cool thermal disk 
might also be detected, with the maximum disk temper¬ 
ature fcTMCD^l keV. The thermal state, which is nor¬ 
mally observed at high luminosities, is characterized by 
a strong disk component with ^Tmcd ~ 0.5-1 keV and 
the disk fractio n of > 7 5 % i n 2-10 keV (as defined in 
[Remillard fc McClintockl 12006 1. A weak Comptonized 
component is also often present but only dominates at 
high energies. The SPL state tends to occur at very 
high luminosities (>0.1 Le<m)> and the hallmark of this 
state is a strong Co mpto n izaton co mpo nent with pho¬ 
ton index Ppl ~ 2.5 dRenrillard fe McClintockl 120061) . A 
sizable thermal component is also normally seen in this 
state. Almost all known Galactic BH LMXBs are tran¬ 
sients. They often show hysteresis in transitions to and 
from the hard state, with the transition to the hard state 
occuring at luminosit ies tig htly c luster ed around 2% of 
the Eddington value dMaccaronel 1200311 . and the transi¬ 
tion from th e hard state at higher but more scattered 
luminosities (iMivamoto et al.l 1 19951 iMaccarone fe Coppil 
120031: iHoman fe Bellonill2005H . Similar hysteresis behav- 
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Fig. 1 . — : Chandra X-ray image of NGC 3115. The im¬ 
age is false-colored using adaptively smoothed (with the 
Cl AO task csmooth), exposure corrected images in 0.5- 
1.2 keV (red), 1.2-2 keV (green), and 2- 7 keV (blue). 
The D 25 ellipse of the galaxy is also shown. The dif¬ 
fuse emission near the galaxy center comes from diffuse 
hot gas, u nresolved LMXBs , and o ther unresolved stellar 
emission (IWong et alJl2nrill2nTl . 

ior is observed in transient NS LMXBs. 

There ar e two mai n classes o f NS LMXBs, at oll and 
Z so urces (lHasinger fe van der Klisl 119891 : Ivan der Klisl 
2006), named after the patterns that they trace out in 
X-ray color-color diagrams (CDs) or hardness-intensity 
diagrams (HIDs). Atolls radiate at ~0.001-0.5 L^dd and 
trace out their patterns in CDs/HIDs on timescales of 
weeks to months. They have two main distinct spectral 
states, i.e., hard (Tpl <2, extending to 100 keV or above) 
and soft states (most emission <20 keV). The hard state 
tends to be observed at low luminosity (<0.1 LEdd)> while 
the soft state is normally observed at high luminosity 
(larger than a few percent LEdd)- Spectra in the “tran¬ 
sitional” state between these two are also occasionally 
seen. Z sources are more luminous than atolls, at near or 
above Eddington luminosity, and they trace out roughly 
Z-shaped tracks in CDs/HIDs within a few days (i.e., 
faster than atolls), with X-ray spectra that are generally 
soft. Thanks to a recent transient Z source XTE J1701- 
462, which exhibited the Z-source characteristics when 
it was accreting at near or above Eddington luminosity 
and transitioned to an atoll during the decay of its 2006- 
2007 outburst, we now know that Z and atoll sources 
are essentially t he same type of objects at different mass 
accretion rates (|Lin et aJ.Jl2009t iHoman et al .11201 Of) . Un¬ 
like BHBs, whose tracks in the CDs/HIDs tend to show 
large scatter, mainly due to the presence of the SPL state, 
NS LMXBs tend to trace out clear narrow tracks in the 


CDs/HIDs, which depend mostly on the accretion rate 
(iDone &; Gierlirisp 2003: R emillajd fe McClintockll2006l : 

Fridri ksson et al.ll2015H . 

The spectral modeling for NS LMXBs is complicated 
by the presence of the boundary layer emissio n and is 
relat ively controversial, compared to BHBs. iLin et akl 
lj2007l) used a similar spectral model to the one used for 
BHBs (described above) except for an additional single¬ 
temperature blackbody (BB) to describe the boundary 
layer and were able to infer L oc T 4 trends for both the 
disk and the boundary layer in the soft state of two atolls. 
Such a trend for the disk is often observed for the thermal 
state of BHBs and is expected if the disk is essentially 
thermal and is trun cated a t the innermost stable circu¬ 
lar orbit (ISCO). ILin et al.l (12009( 1 applied this model to 
XTE J1701-462 and also observed the L oc T 4 trends 
for both the disk and the boundary layer in its atoll 
stage. However, such trends were not observed in its 
Z-source stage, because the inner disk and the boundary 
layer both reach the local Eddington limit in this stage 
so that the increase in the accretion rate tends to lead to 
an increase in the emission area with relatively constant 
temperature (increasing the inner disk radius at constant 
inner disk temperature for the accretion disk). 

Our knowledge of X-ray binaries in nearby galaxies has 
significantly increased since the launch of XMM-Newton 
and Chandra X-ray Observatory in 1999. In particu¬ 
lar, Chand ra’s superb spat i al res olution and excellent 
sensitivity dWeisskonf et al.l [2002H allow for population 
studies of X-ray binaries in a single galaxy with snap¬ 
shots, which is important for understanding the origin 
and evolution of such sources. To fully understand these 
sources, one key task is to identify their nature (i.e., 
BH versus NS X-ray binaries) and the X-ray spectral 
state. To achieve this, early studies tried to stack sources 
in luminosity ranges in order to improv e the statistics 
(|Maccarone et al.l 20031 : llrwin et al.ll2003H . Recently, rel¬ 
atively detailed studies of indi vidual sources in nearby 
galaxies were carried out (e.g.. iBrassineton et al.ll2010L 
120121 : iBurke et al.1120131 : iBarnard et al.lI2013L120141) . but 
limited by statistics, these studies have been mostly fo¬ 
cused on the few most luminous sources (above several 
10 37 erg s _1 ) in each galaxy. It turns out that the dif¬ 
ferentiation between BH and NS X-ray binaries is gen¬ 
erally very difficult. This is because most sources have 
X-ray spectra that are relatively hard (Tpl < 2.5) in the 
Chandra and XMM-Newton bandpass (about 0.3-8 keV) 
and it is difficult to determine whether they are BHBs in 
the hard state or NS LMXBs in the soft state (they are 
unlikely to be NS LMXBs in the hard state due to the 
high luminosities of the sources studied). NS LMXBs in 
the soft state can appear hard in such an energy band 
because of the presence of the h ot boundary layer com - 
ponent in the X-ray spectra (e.g JLm et al.l201(il20l2all . 

NGC 3115 was selected as the target of a 1 Megasec- 
ond Chandra X-ray Visionary Project (XVP) in Cycle 
13. The goals were to study the gas flow inside the 
Bondi radius of the central supermassive black hole and 
obtain a deep look at the X-ray binary population of 
a normal early-type galaxy . The former has been pre¬ 
sented in IWong et al.l (120141) . For the X-ray binaries, we 
have presented the X-ray luminosity function in Lin et 
al. (2014, Paper I hereafter), and here we concentrate 
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TABLE 1 
Observation Log 


Notation 

Obs. ID 

Date 

Exposure 

(ks) 

Offset 3, 

(arcmin) 

1 

2040 

2001-06-14 

35.8 

1.5 

2 

11268 

2010-01-27 

40.6 

0.1 

3 

12095 

2010-01-29 

75.6 

0.1 

4 

13817 

2012-01-18 

171.9 

0.0 

5 

13822 

2012-01-21 

156.6 

0.0 

6 

13819 

2012-01-26 

72.9 

0.0 

7 

13820 

2012-01-31 

184.1 

0.0 

8 

13821 

2012-02-03 

157.9 

0.0 

9 

14383 

2012-04-04 

119.4 

0.3 

10 

14419 

2012-04-05 

46.3 

0.3 

11 

14384 

2012-04-06 

69.7 

0.3 


3 Aim point offset from observation 13820. 


on the detailed properties of discrete sources, such as the 
long-term spectral and flux variability and the spectral 
characteristics. One main goal of our study is to iden¬ 
tify their nature, but different from previous studies, we 
will achieve this by systematic comparison of our sources 
at various luminosity levels with Galactic X-ray binaries. 
NGC 31 15 is a l enticul ar (SO) galaxy with an age of ~8.4 
Gyr dSanchez-Blazouez et, aLll200fij f and at a distance of 
9.7 Mpc (jTonrv et al.ll2001lh Including previous observa¬ 
tions, the total exposure time of Chandra on this galaxy 
is ~1.1 Ms, and the limiting X-ray luminosity is ~10 36 
erg s , making it one of Chandra's best observed normal 
early-type galaxies. 

In Section [2j we describe the source detection, the cal¬ 
culation of flux, the simple spectral fits, the measurement 
of long-term and short-term variability, and multiwave¬ 
length cross-correlation. In Section^ we present the var¬ 
ious properties of X-ray binaries in NGC 3115, including 
the long-term variability and spectral characteristics that 
are used to classify the sources. We further discuss the 
possible nature of our sources in Section [I] We present 
our conclusions in Section [5] 

2. DATA ANALYSIS 

2.1. Observations and Source Detection 

The eleven Chandra observations of NGC 3115 are 
listed in Table [T] They were made during three epochs: 
one in 2001, two in 2010, and nine in 2012. We hereafter 
refer to them as Obs 1-11 in chronological order (Ta¬ 
ble Q]). All observations used the imaging array o f the 
AXA F CCD Imaging Spectrometer (ACIS; iBautz et al.1 
119981 ). We analyzed the data with the Chandra Interac¬ 
tive Analysis of Observations (CIAO, version 4.6) pack¬ 
age. We reprocessed the data to apply the latest calibra¬ 
tion (CALDB 4.5.9) and the subpixel algorithm (|Li et al.l 
I2004D using the CIAO script chandra_repro. Background 
flares are only clearly seen in Obs 1, 5, and 6, for only 
very short durations. We excluded them if they are 
higher than 4cr above the mean background level. In 
this way, we excluded 1.2 ks, 3.6 ks, and 2.7 ks data for 
Obs 1, 5, and 6, respectively. The final exposure used for 
each observation is given in Table [Tj 

To increase the detection sensitivity, we combined all 
11 observations to create a deep merged observation (Obs 
Esum hereafter) using the CIAO script merge-obs. To 
correct for relative astrometry between different obser¬ 
vations, we created new aspect solution files by compar¬ 
ing the source list from each single observation to the 


source list from a single reference observation, which we 
chose to be the longest observation (13820). We only 
used bright sources (>6er) with off-axis angles <6' in the 
cross-correlation. The average separation residuals of the 
source matches are 0.1" after relative astrometry correc¬ 
tion. The new aspect solution files were then applied to 
the event files and subsequent analysis. 

We used the CIAO wavdetect w avelet -based source de¬ 
tection algorithm dFreeman et al.l[2002h to search for dis¬ 
crete X-ray sources. The search was done twice, first 
over the single observations to determine the relative as¬ 
trometry correction described above and the second time 
over the (relative astrometry corrected) single observa¬ 
tions and the merged one. The count images were made 
in the broad (6) energy band 0.5-7.0 keV adopted in 
the Chandra Source Catalog (CSC), while the exposure 
maps were constructed at the corresp onding monochro - 
matic effective energy, (i.e., 2.3 keV. lEvans et al.ll2010 ). 
The point spread function (PSF) maps used correspond 
to the 50% enclosed counts fraction (ECF) at 2.3 keV. 
For the merged observation, the PSF map was obtained 
by averaging those from single observations weighted by 
the exposure. We used two different resolutions: one at 
single sky pixel resolution (0'.'492) over the full field of 
view (FOV) and the other at 1/8 sky pixel resolution 
covering an area of 3'x3' centering at the center of NGC 
3115. The subpixel binning images were used to improve 
the spatial resolution in the crowded field near the cen¬ 
ter of the galaxy. The limiting significance level was set 
to 10 —6 , which formally corresponds to ~1 false source 
due to random fluctuations per 1024x1024 image pixels 
(about one CCD area for images at single sky pixel reso¬ 
lution). The wavelet scales were set to 1, 2, 4, 8, 16, and 
32 image pixels for images at single sky pixel resolution 
and 4, 8, 16, 32, 64, 128, and 256 image pixels for images 
at 1/8 sky pixel resolution. 

Sources detected from the merged observation and 
single observations were cross-correlated to create the 
unique source list. Starting from the source list from 
the merged observation, we searched for new sources 
detected in single observations but not in the merged 
one. Such sources should be faint and highly variable 
or transient. Sources are deemed to be duplicates if 
their separation is less than their combined 99.73% (i.e., 
3cr) statistical positional uncertainty or if their 50% PSF 
circular regions overlap with each other across differ¬ 
ent observations. The statistical positional u ncerta i nties 
that we used are based on Equation (12) of iKim et all 
(|2007f ). which provides the 95% statistical positional un¬ 
certainty as a function of the source net counts and the 
off-axis angle based on a large number of simulations 
using the Chandra simulation tool MARX. The 95% po¬ 
sitional uncertainties are converted to the 99.73% errors 
by multiplying by a factor of 1.405 (this assumed a two- 
dimensional, circularly symmetric Gaussian distribution 
for the source position from wavdetect). As a check on the 
above matching criteria, we also tested a smaller search¬ 
ing radius by using the 95% statistical positional uncer¬ 
tainty only (the 50% PSF overlapping is not used). In the 
end 95 more sources were found. From visual inspection, 
we found that six of them lie in the central crowded field 
and two at the CCD edge, thus probably having large 
systematic positional uncertainties and explaining their 
relatively large offsets from detections from the merged 
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observation. The remaining sources generally have their 
separations from those of the merged observation much 
less than the size of their region ellipse from wavdetect. 
Thus in the end we did not treat any of these 95 sources 
as new sources and used the source list obtained above 
using the 99.73% statistical positional uncertainty and 
the 50% PSF circular region. 

2.2. Flux and Spectral Characterization 

After obtaining a merged source list, we extracted the 
spectrum for each source for each single observation. The 
source region was set to be a circle enclosing 90% of the 
PSF at 2.3 keV. The background region was set to be 
a concentric annulus, with inner and outer radii of two 
and five times the source radius, respectively. Nearby 
sources, if present, were excluded from the source and 
background regions, but the inner circular source region 
enclosing 50% of the PSF was not excluded. We used 
the CIAO task mkacisrmf to create the response matrix 
files and the CIAO tasks mkarf and arfcorr to create 
the point-source aperture corrected auxiliary response 
files. The spectral and response files corresponding to the 
merged observation were created using the CIAO task 
combinespectra. 

The background-subtracted count rates (but not 
aperture-corrected) in different energy bands were ob¬ 
tained from the spectral files. To correctly determine 
confidence bounds for low count limits, we used the 
CIAO task aprates. The conversion from the count 
rates to the fluxes were based on the response files and 
assumed an absorbed PL spectral shape with Tpl = 
1.7 and Galactic absorption Ah = 4.32 x 10 20 cm~ 2 
dKalberla et al.lf2005lh Throughout the paper, all fluxes 
and luminosities quoted (including those obtained from 
spectral fits described below) are corrected for Galactic 
absorption (but not intrinsic absorption, unless indicated 
otherwise). 

To characterize the spectral properties of our sources, 
we calculated the hardness ratios HR = (H — S)/(H+S), 
where S and H are the energy fluxes in the soft and hard 
energy bands, respectiv ely, using the method of Bayesian 
estimation ( Park et all 12006 ). We also carried out sim¬ 
ple spectral fits to spectra above 4cr using two single¬ 
component models: a PL and a MCD. Due to the low 
statistics of most sources, we binned the spectra to have 
a minimum of one count per bin and used the C statistic 
in the fits. Both models included absorption (we used 
the wabs model in XSPEC; we found no significant ef¬ 
fect on our results if we chose other absorption models 
such as tbabs , due to little absorption of most of our 
sources), with the minimum set to be the Galactic value 
of Ah = 4.32 x 10 20 cm -2 . 

2.3. Long-term and Short-term Variability 

The variability of the source was measured in several 
aspects. We defined the long-term flux variability as 
War = Fmax/Fmin and the significance of the difference 
as 

p _ Fmax Fmi n ^, 

- 7^2 +a 2 U/2> Cl 

where F max and Amin are the maximum and minimum 
0.5- 7.0 keV fluxes of a unique source among the single ob¬ 
servations, with the corresponding errors a max and a m i n , 


respectively. We only used detections with the flux above 
twice the error (a) when calculating A max (if no detec¬ 
tions above 2a, the one with the highest significance was 
used as A max ), while we used 2a as the flux for detections 
with the flux less than 2a when calculating A m ; n . 

We measured the shor t-term variability using the 
Gregory-Loredo algorithm dGregorv fc Loredolll992h im - 
plemented by the CIAO tool alvarv (lEvans et al.H 20101. 
It splits the events into multiple time bins and looks for 
significant deviations. The variation of the effective area 
with time was taken into account and was obtained by 
another CIAO tool dither-region. The different degrees 
of confidence are indicated by the parameter of “variabil¬ 
ity index”, which spans values within [0, 1 0] and is larger 
for variability of higher confidence (lEvans et alJl2010f ). 

2.4. Multiwavelength Cross-correlation 

Accompanying the Chandra XVP observation, a six 
pointing Hubble Space Telescope ( HST) mosaic obser¬ 
vation in the F475W and F850LP filters (hereafter g 
and z) using the Advanced Camera for Surveys (ACS) 
was also acquired in the field of NGC 3115. The to¬ 
tal field of view of this mosaic observation is slightly 
larger than the D 25 region of NGC 3115, which has a 
semi-major axis of a = 3.62' (10.2 kpc), a semi-minor 
axis o f b = 1.23' (3.5 k pc) and a position angle of 40° 
lide Vaucouleurs et al.lll99ll ). The galaxy was also im¬ 
aged in g, r, Aband filters on 2008 January 4th using 
Suprime-Cam on the 8.2-m Subaru telescope. In Paper 
I, we have cross-correlated our X-ray sources with the 
360 globular clusters (GCs) from the HST /ACS mosaic 
imaging and the 421 ones from the Subaru/Suprim e-Cam 
imaging (lArnold et al.ll201ll I Jennings et al .1 120 1 III . The 
match was identified if the separation is less than the 
99.73% positional uncertainty (combining the X-ray and 
optical components). The maximum separation allowed 
is 2" in order to limit the spurious rate. We note that be¬ 
fore the cross-correlation, the systematic offset between 
different source lists has been corrected through multi¬ 
ple steps: the Subaru/Suprime-Cam astrometry was reg- 
istered to the USNO-Bl.O Catalog dMonet et al.l 120031 ) 
first, then the HST /ACS astrometry was registered to the 
Subaru/Suprime-Cam one, and in the end the astrometry 
of our X-ray sources was registered to the HST /ACS one 
(therefore the absolute astrometry of our X-ray sources 
has been corrected). 

Tabic [2] lists the 37 matches with HST /ACS GCs (23 
have the g — z color > 1.13 and thus are red/metal-rich, 
while the other 14 have g — z < 1. 13 and are bl ue/m etal- 
poor, following the division in Dennings et all I2014D and 
the 7 matches with Subaru/Suprime-Cam GCs identified 
in Paper I. In Paper I, we also identified five other sources 
wh ose o ptical cou nterparts were not classified as GCs 
by Dennings et alJ (12014D but were assumed to be GC 
candidates by us (Table [2]). Four of them (S12, S53, 
S65, and S79) are within 0.25D 2 5, thus very unlikely to 
be AGNs (Paper I). The other one (S92), at an outer 
region, has an optical counterpart with the size and the 
color typical of GCs, though it has a radial velocity from 
the spectroscopic measurement (238 km s -1 , Table [2]) 
lower than typical values seen in other GCs (> 350 km 

s -1 )- 

Adopting the same matching criteria, we also searched 
for the non-GC counterparts to our X-ray sources outside 
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O. 25 D 25 from these optical observations. Such matches 
are most probably cosmic X-ray background sources 
(CXBs), especially active galactic nuclei (AGNs), in¬ 
stead of the optical counterparts to field LMXBs. The 
optical emission of field LMXBs can achieve the max¬ 
imum when the accretion rate is near the Eddington 
limit and the companion is an evolved star (so bright 
optical emission from both a large disk and a large 
companion) , as seen in the G alactic Z source Cyg X- 
2 (Ivan Paradiis fc McClmtocf] ll995l ). However, such 
sources are still below the detection limit of our optical 
images, by mag. By rotating the X-ray source posi¬ 
tions around the center of the galaxy by ±10°, ±170°, 
and 180° and using X-ray sources above 4er, we estimated 
the rate of spurious matches to be about 3% and 5% for 
the HST /ACS counterparts and Subaru/Suprime-Cam 
matches, respectively. 

From visual inspection, we found that some very faint 
sources in the HST /ACS observ ation were not detecte d 
by the tool SExtractor used by I Jennings et al.l (|2014|) . 
Concentrating on the region within (0.25-l)-D25) we vi¬ 
sually identified five X-ray sources coincident with such 
faint sources and assume them to be AGNs. Within 
0.25-025, we also visually found one so urce, i.e., S65, 
with a faint counterpart not detected bv I Jennings et ahl 
(12014D . and we have assumed it to be a GC, as mentioned 
above. 

3. RESULTS AND DISCUSSION 
3.1. The Source List and Identification 

Figure |T] shows the false-colored Chandra X-ray im¬ 
age of NGC 3115. We detected 525 unique sources from 
the merged and single observations. After eliminating 
sources below 2 a (i.e., the net counts within the 90% PSF 
region divided by the error less than 2), we are left with 
490 sources. We found that the ACIS-S1 chip, which 
is well outside the galaxy (>3_D 2 5), shows some bright 
streaks, especially at energies < 0.7 keV. We eliminated 
seven sources detected in such streaks, as they are most 
probably spurious. We also eliminated another source 
with probably spurious large variability due to being at 
the CCD edge. In the end we have 482 sources in total. 

Table[3]gives the various properties of the sources, such 
as the flux, the long-term variability, and the maximum 
Gregory-Loredo variability index. We also give some 
source type information. For the sources within D 25 , 
they are expected to be dominated by LMXBs and are 
thus assumed to be such objects, except those with non- 
GC optical counterparts, which we identify as AGNs. 
The LMXBs identified in this way are expected to be 
contaminated by some AGNs, because our optical im¬ 
ages are not very deep. Limiting to the 137 relatively 
bright sources with Lx,max > 4 x 10 36 erg s -1 (i.e., flux 
> 3.6 x 10 -16 erg s _1 cm -2 ) within D 25 , which are 
the main targets of our study of this paper, we iden¬ 
tify 9 AGNs. As we found in Paper I that the CXB 
dens ity in our field is consiste nt with the average value 
from lGeorgakakis et al. (l2008t ). to within 20%, we expect 
<16.8 AGNs. Thus the number of AGNs that we might 
miss is < 8 , which is negligibly small for our purposes. 
Sources outside D 25 should be dominated by AGNs, as 
we only expect 1.5 field LMXBs above 4 x 10 36 erg s -1 
based on the IR light in the K s band outside D 25 (Paper 



1 0 37 1 0 38 1 0 39 
kmox (erg s“') 

Fig. 2.— : The long-term variability versus the maximum 
0.5-7 keV luminosity for all LMXBs excluding those in 
the central a = 10" ellipse. The top panel is for GC 
LMXBs, with the blue squares and red diamonds for 
the .ffST/ACS blue/metal-poor and red/metal-rich GCs 
from iJennines et all (12014D . respectively, and the green 
circles for other GCs (the Subaru/Suprime-Cam GCs and 
the ext ra fiv e HST /ACS GC candidates identified by us; 
Section l2~4l) . Two transients, both in red GCs, are plot¬ 
ted with red triangles. The bottom panel is for field 
LMXBs. The blue triangles denote the transient candi¬ 
dates. In both panels, the filled symbols denote BHCs 
(see text for details). 


I). Therefore, we assumed all sources outside D 25 to be 
AGNs, except those coincident with GCs, the supersoft 
X-ray source (SSS) S109, three coronally active stars in 
our Galaxy (they are coincident with stars and show soft 
X-ray spectra and possible stellar flares), two galaxies 
(i.e., due to hot gas emission in galaxies; they are coinci¬ 
dent with galaxies in the optical images and show soft X- 
ray spectra), and the BHC SI79 (Sections 13.21 and 13.3.21 
it is slightly outside D 25 but is also identified as a field 
LMXB, considering its possible transient nature and soft 
spectra). In total we have 136 candidate field LMXBs 
and 49 candiate GC LMXBs (Table [2]). We also marked 
the 13 transients and 10 BH X-ray binary candidates 
(BHCs) in the source type column. Their identification 
will be described in the following sections. 

Table |4] gives the counts, fluxes and hardness ratios of 
our sources in various energy bands in the merged obser¬ 
vation (also in the high-state and low-state observations 
for transients identified in the next section). The counts 
and fluxes of our sources in various energy bands in single 
observations are given in Tabic [5] 
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(a) Transient candidates with the outburst detected in the first epoch. 


n-1... . 


SI 79, V>20 


2074.4 5222 5226 5940 5960 5980 6000 6020 

MJD-50000 (day) 



Energy (keV) 


(b) Transient candidates with the outburst detected in the second epoch. 


Fig. 3.— : Long-term light curve (left panels, with the source number and long-term variability V (i.e., War in 
Section 12751) 1 and sample spectra (right panels, with annotations for the observations and spectral models used) for 
special sources. 


3.2. Long-term Variability 

The observations of NGC 3115 by Chandra span more 
than a decade, which is ideal for investigation of the long¬ 
term variability of LMXBs. Figure [2] shows the depen¬ 
dence of the long-term variability War on the maximum 
0.5-7 keV luminosity Lx,max for all candidate LMXBs ex¬ 
cept those (24) within the central elliptical region with 
semi-major axis a = 10" and eccentricity and position 
angle following the L > 2 5 ellipse (this region is too crowded 
and the source extraction is subject to large systematic 
errors). We find that except for some transients (trian¬ 
gles, see below), War generally decreases with Lx,max for 
both bright GC and held LMXBs with Lx,max > 2 x 10 37 
erg s _1 . At lower luminosities, War seems to increase 


with Lx, max, especially for held LMXBs. However, this 
is most probably artihcial due to the detection limit of 
the observations; most of these sources have the min¬ 
imum luminosity less than 2cr and we calculated their 
War using the 2 a upper limit of the minimum luminos¬ 
ity. Some GCs might contain multiple LMXBs (Paper 
I), but considering that we detect variability for all GC 
LMXBs, such source blending effects should not be sig¬ 
nificant. The long-term stability of the most luminous 
sources that we see in NGC 3115 is also seen in other 
galaxies (e.g.. Ilrwinl 12006 ). 

To search for transients, i.e., quiescent sources with a 
single outburst, we concentrated on the 152 sources that 
are either within 2 D 25 (including those within the cen¬ 
tral a = 10" eclipse) or coincide with GCs and have at 
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(c) Transient candidates with the outburst detected in the third epoch. 

Fig. 3.—: (continue) 
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(d) Transient candidates with the outburst detected in the second and third epochs. 

Fig. 3.—: (continue) 
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(f) A very hard X-ray source (Tpp = 0.6 ±0.1) in a GC, which is a candidate NS LMXB with a strong magnetic 
field or a high inclination. 


Fig. 3.—: (continue) 
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(g) Persistent BH X-ray binary candidates. S96 is in a GC while others are in the field. S108 has only five 
detections because it is in the CCD gap in other observations. 

Fig. 3.—: (continue) 
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least one detection above 4cr. If we required transients 
to have War > 5 (there are 25 such sources), be hardly 
detected (with fluxes < 2 a) in all observations in at least 
one epoch, and have the long-term light curve consis¬ 
tent with a global outburst (instead of an irregular, large 
variation), we are left with 13 candidates. We note that 
we did not use a very strict condition on the variability 
to select transients, as some persistent Ga lactic LMXBs 
are known to vary by factors of >10 (e.g.. iHoman et all 
1 20091 (Maccarone et al.112010 ). Therefore we cannot rule 
out that some transient candidates that we found might 
be just highly variable persistent sources. We plot the 
long-term light curves of the transient candidates in the 
left panels in Figure [3] a-d. There are three active in the 
first epoch (Figure l3a|) . one in the second epoch (Fig¬ 
ure [3b]) , six in the third epoch (Figure [3c]), and thre e ac - 
tive in both the second and the third epoch (Figure |3d|. 
We note that S36 (Figure l3cl). which was bright in the 
third epoch, seems to show some emission in the first two 
epochs, making it not formally a good transient candi¬ 
date, but this might be due to extended emission near 
the galactic center. We also note that only S8 and S103 
are coincident with GCs. 

As most of our observations were made within three 
months in the third epoch, only outbursts in this epoch 
are relatively well monitored. Among the six transient 
candidates that were active only in the third epoch (Fig- 
ure[3c|), three (S103, S125, and S198) clearly showed flux 
evolution during the outburst. We probably detected the 
decay of the outburst for S103 and S125. For S198 we 
had relatively good coverage of the outburst, including 
the rise, the peak and the decay. Moreover, for S198, 
we fortunately caught a fast rise (< 6 d ays), as com¬ 
monly seen in X-ray binaries dChen et al.lll997 ). These 
outbursts over month-long timescales are typical for tran¬ 
sient X-ray binaries. For the other three sources (S8, S26, 
and S36), the flux remained fairly steady in all seven ob¬ 
servations in the third epoch, indicating relatively long 
outbursts (probably years). 

The three transient candidates that were active in 
the second and third epochs had fairly steady fluxes 
in these two epochs. Because the second and third 
epochs together span ~2.2 years, much longer than typi¬ 
cal durations of outbursts seen in X-ray binaries (several 
months), these sources were experiencing prolonged out¬ 
bursts. Such sources are rare but not unique, as there 
are several sources known to ha ve o u tbursts lasting for 
years to more than a decade dMcClintock fc Remillard 
2006: [Galloway et al.l [2007: Sol eri et alJ I2010t iLin et ah 


.2014). 


The sample spectra of all the transient candidates are 
shown in the right panels in Figure [3] Because of large 
variability of these sources and in order to increase the 
statistics, we have combined observations with similar 
flux levels, as noted in the figure. The spectra shown 
generally represent the high state of the sources, except 
S198, for which we also created a low-state spectrum. 
While most spectra appear hard and we show the fits 
with a PL, some spectra seem soft (e.g., S179), and we 
show the fits with a MCD. As will be shown in the next 
section, these soft sources are most probably BHBs in 
the thermal state. We will discuss the spectral evolu¬ 
tion and the physical implication of S198 separately in 
Section 13.4.11 


3.3. X-ray Spectral Properties 
3.3.1. The Hardness Ratios 

Figure [4] shows the luminosity versus hardness ratio 
diagram for our candidate LMXBs, obtained from the 
merged observation (or the merged high-state and low- 
state observations for transients). We separate the GC 
(left panel) and field (right panel) sources. For the for¬ 
mer, we further differentiate different subgroups, i.e., 
the HST /ACS blue/metal-poor GCs in blue squares, the 
HST /ACS red/metal-rich GCs in red diamonds, and the 
Subaru/Suprime-Cam GCs and the extra five HST /ACS 
GC candidates identified by us (Section E3J in green cir¬ 
cles. We observe no clear spectral differences between 
these different groups of GC LM XBs, agreeing with pre¬ 
vious findings (e.g., [Kim et aD 12006 ). For both GC and 
field populations, we find that other than a few very soft 
or very hard outliers, our sources seem to follow a global 
trend that the luminous sources (> 7 x 10 37 erg s _1 ) have 
hard spectra consistent with Tpl ~ 1.5, while the fainter 
sources have systematically softer spectra. In the follow¬ 
ing sections we will present more detailed source spectral 
properties based on simple spectral fits and provide sys¬ 
tematic comparison with Galactic LMXBs, which will 
allow us to shed more light on the nature of our sources 
and the cause of their spectral evolution. 


3.3.2. Spectral fits of field LMXBs 

The results of simple PL and MCD fits to the merged 
spectra (for transients, see Figure [3] for the spectra used) 
of candidate LMXBs are given in TablejS] The C statistic 
that we adopted in the fits does not indicate the fitting 
quality, but based on bright spectra that were rebinned 
to have a minimum of 20 counts per bin and fitted with 
the x 2 statistic, we found that the fits are mostly accept¬ 
able with the reduced y 2 <1.2 and the null hypothesis 
probability > 10%. There are a few relatively bad fits us¬ 
ing the MCD model, with the reduced y 2 of 1.4-1.7 and 
the null hypothesis probability of 10~ 2 10~ 4 . Because 
we are mostly interested in using the PL and MCD fits 
to roughly characterize the spectral shapes, the fitting 
quality overall is sufficient for our purposes. 

The PL fits are shown in Figure El with the 0.3-8 keV 
Lx versus Fpl in the top panels and ./Vh versus Tpl in the 
bottom panels, while the MCD fits are shown in FigureEl 
with Lx versus /cTmcd in the top panels and Ah versus 
^Lmcd in the bottom panels. The GC and field LMXBs 
are plotted separately, with GCs in the left panels and 
field LMXBs in the middle panels. In Figure El some 
dotted reference lines are included to show the expected 
dependence of the 0.3-8 keV Lx luminosity on /cTmcd 
from the standard thermal disk truncated at the ISCO 
of compact objects of several masses (M = 1.4 M 0 , 3 
M 0 , 5 M 0 , 10 M 0 , and 20 M 0 ). We used the empirical 
relation between the real inner disk radi us r; n and the 
normalization of the MCD m odel Amcd ([Kubota et al.l 
H998tlMakishima et al.ll2000h : 


= 1.19A 


' AMCDdi 0k p C 
cos 8 


( 2 ) 


where diokpc is the source distance in units of 10 kpc and 
6 is the inclination angle. The relation takes into account 
the spectral hardening effect, with the hardening factor 
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Fig. 4.— : The 0.5-7 keV luminosity versus the hardness ratio for candidate LMXBs above 4<j (excluding those in the 
central a = 10" ellipse). The hardness ratio is defined as HR = (H — S)/(H + S ), where S and H are the energy 
fluxes in the 0.5-2 keV and 2- 7 keV energy bands, respectively. The luminosities and hardness ratios were calculated 
from the merged observation, except for the 13 transients, for which the merged high-state and low-state observations 
(Figure [3] a-d) were used. The left panel is for GC LMXBs, wi th the blue squares and red diamonds for the HST /ACS 
blue/metal-poor and red/metal-rich GCs from l.Tennings et ahl ( 2014 ). respectively, and the green circles for oth er GCs 
(the Subaru/Suprime-Cam GCs and the extra five HST /ACS GC candidates identified by us; Section l2~4l) . Two 
transients, both in red GCs, are plotted with red triangles. The filled square marks the BHC S96. The right panel is 
for field LMXBs. The blue triangles denote the transient candidates. The filled circles and triangles denote the BHCs. 
The error bars in both panels are at the lcr confidence level. The vertical dotted lines from right to left in both panels 
correspond to hardness ratios for an unabsorbed PL with T =1.5, 2.0, and 2.5, respectively. 


assumed to be 1.7, and the fact that the disk temperature 
does not peak at the inner radius. The dotted reference 
lines in Figure [6] assume 6 = 60°. 

Here we focus on field LMXBs first, with GC LMXBs 
to be presented in the next section. The most striking 
result of the PL fits to field LMXBs is the strong depen¬ 
dence of the photon index on the luminosity, as shown in 
the top middle panel in Figure [5] Most sources fall on a 
narrow track (the light gray region), with Tpl decreasing 
(thus the sources becoming harder) with increasing Lx 
up to Lx ~ 7 x 10 37 erg s _1 and then remaining at a 
value around 1.5 above this luminosity. Such a trend has 
been indicated, though with larger scatters, in the lumi¬ 
nosity versus hardness diagram in Figure 0] A few very 
soft outliers lying on the right of the light gray track can 
also be seen. Five of them are persistent sources (filled 
circles), and the other four are transients (filled triangles; 
S198 has two data points, corresponding to its high and 
low states). At Lx < 2x 10 37 erg s _1 , some hard sources 
with r pl < 1.8 are also present. We note that the lack 
of hard sources with Fpl < 1.8 at Lx < 3 x 10 36 erg s _1 
is due to selection bias because we only fitted sources 
above 4cr and harder sources tend to have lower signifi¬ 
cance levels at a given luminosity. The column density 
was inferred to be < 0.4 x 10 22 cm -2 for most sources. 
It seems to increase with the photon index, which could 
be caused by the use of a (non-physical) PL model to 
fit the spectra that are mainly thermal and are softer at 
lower luminosity (more discussion on the column density 
is given at the end of the section). 

The MCD fits give results consistent with the PL fits. 


Most sources reside in the light gray region below the 
M = 3 Mq dotted reference line (thus corresponding to 
lower mass) in the top middle panel plotting Lx verus 
kTucn in Figure 0 with kTyiGD increasing with Lx be¬ 
low ~7 x 10 37 erg s -1 and then remaining at around 1.5 
keV above this luminosity. The very soft outliers identi¬ 
fied from the PL fits above fall in the region correspond¬ 
ing to M > 5 M 0 . The column density of most sources 
from the MCD fits is at the minimum value allowed in 
the fits, i.e., the Galactic value. 

To shed light on the possible nature of the sources, in 
the right panels in both Figures 0 and 0 we plot the ex¬ 
pected single-component (PL and MCD) fitting results 
of three representative Galactic X-ray binaries, the atoll 
source 4U 1705-44 (filled red diamonds), the Z source 
GX 17+2 (green crosses) and the BHC XTE J1817- 
330 (filled blue squares), based on the spectral fi ts of 
these sources bvlLin et al.1 (|2010H . lLin et ahl (l2012alh and 
iRvkoff et alJ (|2007f) . respectively, and assuming them to 
be at the distance of NGC 3115 with absorption at the 
Galactic value in the direction of NGC 3115 (i.e., as¬ 
suming zero intrinsic absorption). One main reason for 
us to choose these sources to compare with our sample 
is that they had high-quality broad-band spectra span¬ 
ning large luminosity ranges and were carefully mod¬ 
eled in the ab ove studies. The spectra of 4U 1705-44 
studied by iLin et all 1 2010 1 used broad-band observa¬ 
tions by Suzaku (1.2-40 keV) and BeppoSAX (1 -150 
keV) and included two hard-state spectra and seven soft- 
state spectra over a l arge dynamical ra nge (~0.04-0.28 
LEdd)- For GX 17+2. ILin et al.l (l2012aT) compiled spec- 
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Fig. 5.— : The PL fit results of LMXBs excluding those in the central a = 10" ellipse, with the top panels plotting 
the 0.3-8 keV luminosity versus the photon index (symbol size proportional to column density logarithm) and the 
bottom panels plotting the column density (including the Galactic absorption) versus the photon index (symbol size 
proportional to luminosity logarithm, using data with Lx > 8 x 10 36 erg s -1 for clarity). The left and middle panels 
are for GC and field LMXBs, respectively (the meanings of the symbols are the same as in Figure |4}. The right panels 
plot the PL fits to the atoll source 4U 1705—44 (filled diamonds), the Z s ource GX 17+2 (crosses) a nd the BHC 
X TE J181 7—330 (filled squares), based on the spectral fits of these sources bv lLin et al.l (12010T) . ILin et all (l2012alh and 
iRvkoff e~a l. (2007), respectively, and assuming them to be at the distance of NGC 3115 with only Galactic absorption 
(jV h = 4.32 x 10 20 cm -2 ). In these panels, we also plot the data shown in the middle panels for field LMXBs, but in 
a gray color. The error bars (sometimes smaller than the symbol size) in all panels correspond to the 90% confidence 
level. The light gray region in the top panels marks the possible NS LMXB soft state track, where atolls in the soft 
state and Z sources reside. 


tra over the whole Z track in the HID and used RXTE 
data (~ 3-40 keV), but for clarity we only used spec¬ 
tra at four critical positions: the two vertices and the 
two ends of the Z track in the H ID . W e select XTE 
J1817-330 studied bv IRvkoff et al.l (|2007H to represent 
the spectral properties of BHBs because the outburst 
in their study was covered by Swift/XKT in the soft 
X-ray energy band of 0.5-10 keV, close to the Chan¬ 
dra bandpass, and spanned a large dynamical range, 
while the source is also one of the Galactic BHCs with 
the lowest absorption (IVh = 6 x 10 2 ° cm -2 ). Several 
spectral models were applied to these three sources in 


the above three studies, but we used the fits of the 
models that carefully addressed Comptonization, i.e., 
the SIMPL(MCD)+BB model, the MCD+B B+n thcomp 
model and the M CD+ comptt model in ILin et all (120101) . 
ILin et all (l2012af) . and IRvkoff et al.l (120071) . respectively. 
This is important because Comptonization with a steep 
slope can significantly overestimate the low-energy emis¬ 
sion if it is fitted with models like a simple PL. 

To convert the source flux to that at the distance 
of NGC 3115, we have ado pted the source d ist auce o f 
d = 7.4 kpc for 4U 1705-4 4 (lHaberl fc Titarchuklil995 ). 
d = 8.1 kpc for GX 17+2 (ILin et al.ll2012afl . and d = 10 
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0.1 1 0.1 1 0.1 1 

kT MCD (keV) kT MCD (keV) kT MCD (keV) 

Fig. 6.—: Similar to Figure[5]but for the MCD fits. The dotted lines in the upper panels show the expected dependence 
of the 0.3-8 keV luminosity Lx on kTyicu from the standard thermal disk truncated at the ISCO of compact objects 
of several masses (M = 1.4 M 0 , 3 M 0 , 5 M 0 , 10 M 0 , and 20 M 0 corresponding to lines from the right to the left, 
respectively; see text for details. The light gray region in the top panels marks the possible NS LMXB soft state track, 
where atolls in the soft state and Z sources reside. 


kpc for XTE J1817-330 (as assumed in [Rvkoff et al.l 
12001 . The distance of 4U 1705-44 was derived from 
type I X-ray bursts and has a relatively small uncertainty 
(~20%) dependi ng on the abundance of the accreted ma¬ 
terial assumed (|Gallowav et al.l [20081) . The distance of 
GX 17+2 was also derived based on type I X-ray bursts, 
but the bursts in this source exhibited abnormal prop¬ 
erties, which could lead to a rel atively large uncertainty 
(~50%) on the distance derived dLin et al.ll2012af) . There 
is no reliable distance estimate for XTE J1817-330. The 
distance of 10 kpc impli ed a BH of 10 M 0 , based on 
Equation [2] dRvkoff et al.l 1200 7l) . As long as the central 
BH is not too small and the Eddington ratio of the ther¬ 
mal state in this source is not too low (> 1%), we expect 
its distance to be close to 10 kpc (within a factor of 2). 
The above distance uncertainties of GX 17+2 and XTE 
J1817-330 are small enough and do not significantly af¬ 
fect our comparison of these sources with the sources in 
NGC 3115. 


For each spectrum of the representative Galactic X-ray 
binaries, we generated 200 simulations and fitted them 
with PL and MCD models. The data plotted in Figures[5] 
and [S] are the median of the fitting results, with the error 
bars representing the intervals including 90% of the fits. 
In the PL and MCD fits to the simulated 0.3-8 keV spec¬ 
tra of the seven soft-state observations of 4U 1705-44, we 
inferred the photon index to decrease from 1.9 to 1.3 (the 
top right panel in Figure [5]) or the disk temperature in¬ 
creasing from 1.1 keV to 2.1 keV (the top right panel in 
Figure[6|) as the 0.3-8 keV luminosity increases from 10 37 
erg s _1 to 7 x 10 37 erg s _1 . Such trends are very similar 
to those traced out by a majority of our field LMXBs 
in the light gray region in Figures [5] and [6] in a similar 
luminosity range. GX 17+2 shows similarly hard spectra 
to the brightest observation of 4U 1705-44* in 0.3-8 keV, 
with Fpl ~ 1.4 from the PL fits and &Tmcd ~ 2 keV 
from the MCD fits. Such properties are similar to those 
of the brightest held sources. 
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Based on the above comparison, we suggest that the 
light gray region in Figures [5] and E] marks a NS LMXB 
soft state track, with sources below Lx ~ 7 x 10 37 erg s _1 
being atolls in the soft state and brighter sources be¬ 
ing Z sources. Then our observation that most of our 
sources appear harder at higher luminosity in 0.3-8 keV 
can be easily explained because for the atoll soft state, 
the temperatures of the disk and boundary layer ther¬ 
mal emission increas e with luminosity at relative ly con¬ 
stant emission areas (IT, in et al .1120071 [20M [2tTmh . Such 
trends are not seen in Z sources because the inner disk 
and the boundary layer reach the local Eddington limit, 
in which case the change in the accretion rate tends to 
lead to change in the emission area (the in ner rad ius for 
the disk), instead of the temperature dLin et al.l I2009H . 
We note that our sources seem systematically slightly 
softer than the simulated spectra of Galactic NS LMXBs, 
and one explanation for this is the possible presence of 
soft excess of the real disk spectrum compared with the 
si mple MCD desc riptio n (i.e., diskbb in XSPEC) used 
in ILin et ahl (l2010f) and iLin et al~l (|2012af) , which will be 
discussed in Section FOl 

Although the soft-state observations of 4U 1705-44 in 
ILin et al~ (j2010T) are above Lx ~ 10 37 erg s _1 , atolls 
in the soft state can be fainter. For example, the soft 
state of 4U 1608-52 has Eddington ratios as low as 0.01 
(ILin et al.ll2007lh four times lower than 4U 1705-44 in 
ILin et aLpl2010f ). Therefore, the atoll soft-state explana¬ 
tion could apply to our relatively soft (Tpl > 2.0) faint 
sources with Lx as low as a few 10 36 erg s _1 . 

As mentioned above, we also have nine very soft out¬ 
liers (filled circles and triangles in the middle panels in 
Figures [5] and 0- They are much softer than NS LMXBs 
at corresponding luminosities, but are similar to the BHC 
XTE J1817-330 in the thermal state. Therefore, we iden¬ 
tified them as BHCs. Their long-term luminosity curves 
and sample spectra are shown in Figure [3] (those fitted 
with a MCD or a MCD+PL model), and the PL and 
MCD fit results are given in Table [6] The spectra of 
most of these BHCs can be described with a MCD, ex¬ 
cept S100. This source seems to show a hard tail, and 
we tried to fit it with a MCD plus a PL, with the pho¬ 
ton index fixed at a value of 2.5. The fitting result is 
included in Table [71 indicating the presence of a PL at 
the 6cr confidence level. S198 is an interesting BHC that 
will be described separately in Section 13.4.11 It is the 
only source with two spectra (one in the high state and 
the other in the low state, Section 13.4.11 Figure [3c]) in 
Figures [5] and El 

The faint hard sources with Lx < 2 x 10 37 erg s _1 and 
Fpl ^ 1.8 could be atolls in the hard state, similar to 
the two hard-state spectra of 4U 1705-44 (the two filled 
diamonds below 8 x 10 36 erg s _1 in Figures E] and E]). 
We cannot rule out the possibility that there may be 
some BHBs in the hard state. For example, simulating 
the h ard-state s pectr um of the BHB GX 339-4 studied 
by iMiller et al.1 (1200611 and fitting with a PL, as we did 
above for other three Galactic X-ray binaries, gave Lx ~ 
2 x 10 37 erg s 1 and Tpl ~ 1.6. However, considering 
that NS LMXBs are expected to be more common than 
BHBs, the possibility being BHBs should be low. 

Our PL fits to simulated spectra of three Galactic X- 
ray binaries in the soft/thermal state systematically gave 


column densities larger than the value assumed in the 
simulation (i.e., the Galactic value toward NGC 3115; see 
the bottom right panel in Figure [5]). Therefore, the col¬ 
umn densities from the PL fits to our sources are system¬ 
atically overestimated if our identification of our sources 
(i.e., most sources are NS LMXBs in the soft state with 
the very soft outliers being BHBs in the thermal state) is 
correct. On the contrary, the MCD fits tend to underes¬ 
timate the column density, due to the presence of a hard 
component (weak Conrptonization in the case of BHBs in 
the thermal state or the boundary layer emission in the 
case of NS LMXBs in the soft state) and/or possible ex¬ 
tra soft emission of the real disk spectrum relative to the 
MCD description, which will be discussed in Section [4] 

There is one main caveat of the above comparison of 
our sources with Galactic sources to be kept in mind. 
To have enough statistics, we have to use the spectra 
of our sources accumulated over long exposures (~1.1 
Ms) and spanning more than a decade, though most 
data were made from the XVP, which still covered nearly 
three months. To reduce this problem, we have created 
spectra for transients carefully by combining only ob¬ 
servations with similar fluxes. For the other sources, if 
Lx, max > 2 x 10 37 erg s _1 , most sources have the long¬ 
term variability V val < 3, which is not large enough to 
affect our comparison with Galactic sources significantly. 
However, for fainter sources, their long-term variability 
is uncertain, and there might be relatively large system¬ 
atic errors produced by the dependence of the spectral 
properties on the luminosity obtained. 

3.3.3. Spectral fits of GC LMXBs 

On the whole, the simple PL and MCD fits indicate 
that GC LMXBs have spectral properties very similar 
to those of field LMXBs (compare the left panels for GC 
LMXBs and the right panels for field LMXBs in Figures[5] 
and El). In particular, most GC sources are also in the 
NS LMXB soft state track (the light gray region), which 
exhibits a trend of decreasing Tpl or increasing kTyicu 
with increase in Lx below Lx ~ 7 x 10 37 erg s _1 and 
relatively constant Tpl (~ 1.5) or AiImcd (~ 1.5 keV) 
above this luminosity. As discussed in the previous sec¬ 
tion for field LMXBs, most of our GC LMXBs below and 
above Lx ~ 7 x 10 37 erg s _1 in this track could be atolls 
in the soft state and Z sources, respectively. Some hard 
sources with P < 1.8 below Lx ~ 3 x 10 3 ' erg s _1 are also 
present, and they could be NS or BH LMXBs in the hard 
state. It seems that GC LMXBs have more hard sources 
at luminosities around 10 37 erg s _1 than held LMXBs. 
For example, over the range Lx = (1-2) x 10 37 erg s _1 , 
there are six out of 10 GC LMXBs with T < 1.6, but 
there is only one out of 14 held LMXBs with T < 1.6, 
from the PL hts. 

There is one persistent source S96 in a blue/metal-poor 
GC (the filled square in the left panels in Figures Eland© 
whose spectrum is signihcantly softer than those of other 
sources at similar luminosity (Lx ~ 3 x 10 37 erg s" 1 ). 
Its long-term luminosity curve and merged spectrum are 
shown in Figure [3gJ We identihed it as a BHC. There 
is another source S223 much harder than other sources 
at similar luminosity (Lx ~ 10 38 erg s _1 ), which is a 
candidate NS LMXB with a high magnetic held or a 
high inclination and will be presented separately in Sec¬ 
tion [3X31 
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3.3.4. CXB sources 

The results of PL fits to CXB sources are shown in 
Figure [7] Sources inside and outside 2D 2 5 are plotted in 
the left and right panels, respectively. For both groups of 
sources, we see no clear dependence of the photon index 
on the flux and column density. Using sources outside 
2Z?25 and excluding extreme sources with Tpl > 3.5 and 
those with Tpl < 0.5 and limiting to Lx > 2 x 10~ 15 
erg s -1 cm -2 , we obtained a median of Tpl = 1.77 with 
68.3% within 0.3 8 ar ound it, consistent with the results 
from lLin et all (|2012bT> . 

In Figure H we include three sources (triangles; S317, 
S362, and S372) that are candidate coronally active stars 
because they all have Tpl > 3.5, coincide with bright 
stars in the Subaru/Suprime-Cam imaging, and have 
possible stellar flares detected, as shown in the bot¬ 
tom three panels in Figure [8] Two faint soft sources 
(T pl > 2.6; S288 and S290) are plotted with squares. 
Since they are coincident with bright extended galaxies 
in the Subaru/Suprime-Cam imaging, they might be due 
to hot gas emission in the galaxies. 

3.4. Special Sources 

3.4.1. SI98: A transient BH LMXB candidate with a 
prominent disk at a very low state 

As shown in Section [T2l S198 is the only source whose 
outburst is covered relatively well by our observations 
f Figure l3cl). It also has the most clear spectral evolu¬ 
tion caught among our sources. To accumulate enough 
statistics, we combined all the data in the rising (ob¬ 
servations 4-5) and decay (observations 9-11) stages to 
create a low-state spectrum (the spectra in these two 
stages seem to be consistently soft) and the others (6- 
8) for a high-state spectrum. When we adopted the PL 
model, we obtained Tpl = 2.6 ± 0.3 for the high state 
and Tpl > 4.2 (the 90% lower limit) for the low state 
(the absorption for the low state was fixed at the value 
obtained from the fit to the high state; Table HJ). Us¬ 
ing the MCD model, we obtained /cTmcd = 0.78 ± 0.07 
keV for the high state and /cTmcd = 0.17 ± 0.05 keV for 
the low state. Therefore, the spectrum of S198 in the 
low state in the 2012 outburst is much softer than in the 
high state. The 0.3-8 keV Lx is about 0.12 and 0.008 
Eddington luminosity (assuming a BH of mass Mbh=6 
Mq (see below) and a disk inclination of 60°) in the high 
and low states, respectively. 

The disk temperature and luminosity in the high 
state are typical for BH X-ray binaries in t he t hermal 
state (Figure [6] also refer to Figure 16 in iDone et al.1 
f200l . in which the thermal _ disk emission dominates 
i Remillard fc McClintockl [2006 '). Based on Equation [2] 
and the MCD fit in Table [6] and assuming a disk inclina¬ 
tion of 60°, the BH has a mass about 6 Mq. 

The low state is more difficult to understand, and 
to determine the corresponding spectral st ate, we com ¬ 
pare our source with Swift J1753.5-0127 in iRvkoff et al.1 
( 2007 ?). We first checked for the possible presence of a PL 
component in the low state by adding a PL to the MCD 
model in the fit to the spectrum. We fixed Tpl = 1-7 and 
Tpl = 2.5, typical values seen in BHBs, and found that 
the PL component contributes less than 23% and 37% 
(90% upper limit) of the 0.3-8 keV unabsorbed flux, re¬ 
spectively. Thus the MCD still dominates in this energy 


range in such fits. In comparison, Swift J1753.5-0127 has 
a disk thermal flux fraction of about 71%, 62%, and 56% 
in observations 15, 16, and 17 (those with Lx ~ 1.4xl0 3 ' 
erg s _1 , 2.5 x 10 3 7 erg s -1 , and 3 .4 x 10 37 erg s _1 in Fig¬ 
ures [5] and [6| in IRvkoff et al.l (1200711 , respectively. The 
disk temperatures are about 0.19 keV, 0.25 keV, and 0.33 
keV, and the Eddington ratios are about 0.010, 0.019, 
and 0.026, assuming a BH mass of 10 Mq and a source 
distance of 10 kpc, in these three observations, respec¬ 
tively (IRvkoff et al.l l2007il . The source was probably in 
the transitional state in observations 15 and 16 and in 
the hard state in observation 17, based on the broad¬ 
band data from the Proportional Counter Array onboar d 
the R ossi X-ray Timing E xplorer (IGierlinski et all2008t ). 
While IRvkoff et all (|2007f ) concluded that the inner disk 
radius was consistent with the IS CO in all the three and 
other brighter Swift observations. IGierlinski et al.1 ([20081 1 
demonstrated that the inner disk receded from the ISCO 
in observations 15-17, after taking into account the effect 
of irradiation from the hot corona on the disk. The inner 
disk radius in the low state of S198 is about three times 
larger than that in the high/thermal state, though only 
at the 2.2cr level (Table [6]). Therefore, based on the disk 
temperature, the Eddington ratio and the larger inner 
disk radius than in the high state, S198 in the low state 
could be similar to Swift J1753.5-0127 in observation 17, 
i.e., in the hard state. Alternatively and probably more 
likely, it could be in the transitional state, considering 
its more prominent disk than Swift J1753.5-0127 in ob¬ 
servation 17. 

3.4.2. S109: A SSS 

SSSs have a characteristic temperature <0.1 keV 
and have been seen in our Milky Way, the Magellanic 
Clouds, and nearby g alaxies (for a recent review, see 
iDi Stefano et al.1 120101) . S109 is the only SSS in our 

source sample. Its long-term light curve and sum spec¬ 
trum are shown in Figure [3(3 The source seems persis¬ 
tent with lUar = 3. The spectrum can be fitted with 
a BB with /cTbb = 86j^ 2 eV, apparent emission radius 
R B b = 3.41q 3 x 10 3 km (TablcjTJ). The unabsorbed bolo- 
metric luminosity is (O.S+g^) x 10 38 erg s _1 . SSSs with 
luminosity around the Eddington limit for a white dwarf 
(WD), i.e., ~ 10 38 erg s _1 , coul d be n uclear burning of 
material accreted by a WD (lGreinedl2000f l. Therefore 
our SSS is probably one such source. 

3.4.3. S223: A very hard luminous X-ray source in a GC 

As shown in Section 13.3.31 S223 has a spectrum sig¬ 
nificantly harder (Tpl = 0.6 ± 0.1) than other sources 
(r pl > 1.4). The source seems persistent, with a long¬ 
term variability factor of V var = 1.8 at the 2.2a signifi¬ 
cance level and a mean 0.3-8 keV luminosity of 0.96x 10 38 
erg s _1 (Figure (33) ■ We detected short-term variability 
in observation 11, with a Gregory-Loredo variability in¬ 
dex of 6 (Figure [Sj, but not in other observations, with 
Gregory-Loredo variability index of < 1. The source is 
coincident with a spectroscopically confirmed GC, whose 
g — z color (1.1) is close to the boundary value of 1.13 
used to separate the blue and red GCs (Paper I). There¬ 
fore it is very unlikely to be a background AGN or a 
foreground star. It is probably not composed of multi¬ 
ple sources in the same GC, considering the detection of 
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Fig. 7.—: Similar to Figure [5] but for the PL fits to candidate CXB sources. Sources inside and outside 2D 2 b are 
plotted in the left and right panels, respectively. The bottom panels for Nn versus Tpl only show sources with 
Lx > 8 x 10 -16 erg s -1 cm" 2 . Three candidate coronally active stars (triangles in the right panels) are also included. 
Two faint soft sources (squares in the right panels) could be from hot gas emission of optically bright galaxies. 


both long-term and short-term variability. 

There are other X-ray sources that were found to show 
very hard spectra and high luminosities and reside in 
GCs in nearby galaxies. IBurke et al.l (I2013T) reported 
a source (their S48) in a GC in Centaurus A having 
Tpl ~ 0.7 and the 0.5- 10 keV luminosity ~ 6 x 10 37 erg 
s _1 . iTrudolvubov fe Priedhorskvl (120041 ) reported two 
sources (their S22 and S32) in two GCs (Bo91 and Bol58, 
respectively) in M31 having Tpl ~ 0.6-0.8 and the 0.3- 
10 keV luminosity ~(0.6-1.8)x 10 38 erg s _1 . These two 
sources showed a spectral cutoff of ~4-8 keV (our source 
showed no cutoff below 12 keV at the 90% confidence 
level). It is not_ cle ar whe ther these sources have a com¬ 
mon nature. I Burke et~ . (1201 34 suggested their S48 as a 
highly magnetized NS as seen in normal accreting X-ray 
pulsars (instead of millisecond X-ray pulsars), which typ¬ 
ically have r PL < 1, cut-off e nergies ar ound 20 keV, and 
luminosities > 10 36 erg s _1 (White et al.l 11983ft . Con¬ 
sidering that such objects normally show coherent pul¬ 
sations, we searched for them for our source by creat¬ 


ing power densi ty spectra for each observation (refer to 
iLin et all (2014 ) for the procedure). We found no detec¬ 
tion above the 99% significance level in any observation 
from the timescale of the whole observation length up 
to the Nyquist frequencies (0.62 Hz for observations 1-3 
and 0.64 Hz for observations 4-11). It is possible that its 
pulsation period is shorter than the readout frame time 
or that the pulsation is not strong enough to be detected 
with current data. Most normal accreting X-ray pulsars 
are high-ma ss X-ray binaries, with only a few known to 
be LMXBs (iBildsten et al.lfT997D . Given the coincidence 
with a GC, S223 is more likely to be a LMXB if it is a 
normal accreting X-ray pulsar. 

We note that the hard X-ray source in Bol58 in 
M31 showed periodic dips oc casionally and thus prob- 
ably ha s a high inclination (ITrudolvubov et ahl 120021 
ITrudolvubov fc Priedhorskvll2004f ). Therefore we specu¬ 
late an alternative explanation for the hard spectra of the 
above sources: they might have high inclination angles, 
which suppress the observed disk emission and enhance 
the observed boundary layer emission and thus cause the 
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Fig. 8.— : The light curves of some sou rces with in¬ 
teresting variability, with Gehrels errors (lGehreldll986il 
shown. The top panel is for S223, which is a candi¬ 
date NS LMXB with a strong magnetic field or a high 
inclination in a GC, in observation 11, showing signifi¬ 
cant short-term variability. The bottom three panels are 
for three coronally active star candidates in observations 
showing sign of flares. The bin size At of each light curve 
is annotated in each panel. 

hard X-ray spectra (but the inclination should not be too 
high for the boundary layer to be strongly obscured by 
the disk). Limited by statistics, we cannot determine 
whether the short-term variability of our source S223 in 
observation 11 is due to dipping. 

4. DISCUSSION 

4.1. The NS LMXB Soft State Track and The New 
Source Identification Scheme 

Except some very soft outliers, which should be strong 
BHCs in the thermal state, the majority of our bright 
sources appear hard in the Chandra bandpass (0.3-8 
keV). They could be BHBs in the hard state or NS 


LMXBs in the soft state. The spectra of NS LMXBs 
in the soft state can appear hard in this energy band 
due to the hot boundary layer emission. Differentiating 
between the above two scenarios is nontrivial based on 
the narrow-band spectra in hand. The method that we 
adopted is to compare the collective spectral properties 
of our sources with the spectral evolution of three repre¬ 
sentative Galactic X-ray binaries based on simple PL and 
MCD fits. We found that most of our sources fall on a 
narrow track in the Lx versus Tpl and Lx versus fcXivrCD 
plots, exhibiting harder spectra at higher luminosity be¬ 
low Lx ~ 7 x 10 37 erg s -1 but relatively constant spectral 
shape above this luminosity. Such spectral evolution is 
close to that expected for NS LMXBs in the soft state in 
the Chandra bandpass. Therefore, we identify the track 
as the NS LMXB soft state track, in which sources below 
Lx ~ 7 x 10 37 erg s _1 are most likely atolls in the soft 
state and sources above this luminosity are Z sources. 

Our sample of candidate LMXBs in NGC 3115 also in¬ 
cludes some hard sources at low luminosities. Although 
we believe that they should be dominated by NS LMXBs, 
the possibility of some being BHBs (even AGNs) cannot 
be ruled out. Therefore our list of BHBs should be con¬ 
servative, only including the ten strong BHCs that were 
identified based on their very soft spectra. Restricting 
our search to (0.046-1.0)D 2 5 and above 10 37 erg s _1 , we 
found one persistent BHC out of 23 GC LMXBs and 
nine BHCs (five persistent and four transient) out of 
59 field LMXBs. Therefore, we found a much larger 
fraction of BHBs in the held than in GCs, as seen in 
our Galaxy. We have identified four BHCs out of 11 
held tra nsients . This is more abundant than the simula¬ 
tions bv lFraeos et al.l (12008D (<10% for transients above 
10 37 erg s _1 ). Our discovery of hve persistent BHCs in 
the held with Lx > 10 37 erg is also not predicted by 
iFragos et al.l ( 2008 ). 

X-ray sources in the old populations in other galax¬ 
ies seem to show similar spectral properties to those of 
our sources a nd can thus be i dent ified in the same way. 
iTrudolvubov fc Priedhorskvl (j2004f l htted the XMM- 
Newton (using the 0.3 10 keV band) and Chandra (using 
the 0.5-7 keV band) spectra of 31 bright GC LMXBs in 
M31 with an absorbed PL. Their sources had the 0.3- 
10 keV luminosity Lx in the range of ~ 10 36 erg s _1 to 
10 39 erg s _1 . In the Lx versus Tpl plot (their Figure 
6), we can see that most of their sources reside in the 
NS LMXB soft state track and thus should be atolls in 
the soft state and Z sources following our identification 
scheme, except some hard sources at low luminosities, 
which could be NS (most li kely) or BH LMXBs in the 
hard state. ITrudolvubov fe Priedhorskvl (j2004f) also ar¬ 
gued that most of their sources, almost all being per¬ 
sistent, should be NS LMXBs, based on the similarity 
of their spectral properties and long-term variability to 
those of luminous persistent X-ray binaries in our Galaxy 
(they are mostly NS LMXBs). However, in their fits to 
the three brightest (Lx > 10 38 erg s _1 ) sources with 
the MCD+BB model, two have /cTmcd 1$ 1 keV, which 
seems too low for Z sources, as will be discussed below. 

In the stud y of X-ray bin aries i n Centaurus A observed 
by Chandra , iBurke et ali (i2013h also obtained the Lx 
versus Tpl plot and the Lx versus AiTmcd plot (their 
Figures 3 and 2, respectively), but for only some bright 
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sources (0.5-10 keV luminosity Lx in the range of ~ 10 37 
erg s _1 to 4 x 10 38 erg s _1 ). They tried to determine 
whether the source spectra were dominated by a MCD or 
a PL based on the behavior o f the inferred column den - 
sity. Their method followed iBrassington et all (|201Q[ ). 
who came up with this method based on simulations of 
typical BHB spectra. We prefer our source identification 
scheme, which is based on the spectral behavior of both 
NS and BH LMXB s. Almost all the 1 7 sources in the Lx 
versus Tpl plot in lBurke et al.1 (120131 1 that they found to 
be PL dominated are in the NS LMXB soft state track. 
Among the 1 6 sou rces in the Lx versus /cTmcd plot in 
IBurke et al.1 (|2013( 1 that they found to be MCD domi¬ 
nated, 12 are in the NS LMXB soft state track, while 
the other 4 have very soft spectra (/cTmcd ^_1 keV) and 
are pr obably BHBs in the thermal state. IBurke et al.1 
(12013!) generally gave consistent source identifications, 
mostly based on the disk temperature from the MCD fits 
only, and our more detailed comparison with Galactic NS 
LMXBs makes the identification of the source nature and 
spectral state more convincing. 

4.2. Comparison with The Double Thermal Model 
Source Identification Method 

Given that ILin et al.l (120071 120091 l20l(il I20l2all suc¬ 
cessfully fitted thousands of NS LMXB soft-state spectra 
with the double thermal model (MCD+BB), pl us a weak 
Comptonized component when necessary, iBarnard et al.1 
( 2013 1 applied the double thermal model to 35 sources 
in M31 showing bright hard spectra and found them to 
be BHCs because their best-fitting parameters lie sig¬ 
nificantly outside the space occupied by NS LMXBs. In 
particular, they found that the disk temperatures of their 
BHCs from the fits with the double thermal model are 
systematically lower (/cTmcd <1.0 keV) than those of 
typical NS LMXB soft-state spectra (/cTmcd <1 1.0 keV), 
and the BB fractions of the 2-10 keV f l ux are also larger 
(>65% versus <50%). IBarnard et al.l (I2014T) expanded 
the work and found 15 more BHCs, i.e., 50 in total. 
These 50 B HCs incl ude 15 GC so u rces , 12 of which were 
included in iTrudolvubov fc Priedhorskvl (20041 1. Based 
on our source identification scheme, however, these 15 
BHCs (two are transients) in GCs should be NS LMXBs. 
For the other 35 non-GC sources, most of them should be 
NS LMXBs too, except for about six sources with very 
soft spectra, which should be BHCs. 

We also tried to fit the double thermal model to our 
bright sources with Lx > 4 x 10 3 ' erg s _1 . There are 
33 (14 in GCs and 19 in the field), excluding the soft 
BHCs that we have identified. We found that they all 
have /cTmcd < 1 keV, and 27 of them have the upper 
error o f /cTmcd at the 90% c onfid ence level < 1 keV. Fo l- 
lowing IBarnard et al.l (|2013l ) and IBarnard et al.1 (|2014H . 
this would indicate that almost all bright sources with 
L x > 4 x 10 37 erg s” 1 in NGC 3115 are BHCs, which 
seems very unlikely, considering that BHBs are expected 
to be much rarer tha n NS LM XBs in an early-type galaxy 
like NGC 3115 dFragos et al. 1120081) . 

We note that we are using a narrower and softer energy 
band and have sources sub j ect to much less absorption 
than ILin et al.1 (120071 120091 1201(1 l2012alf . These differ¬ 
ences could cause problems if the double thermal model 
is directly applied to our sources. We estimated the pos¬ 
sible systematic errors arising from the use of a narrow 


<N 

O 

O 


CM 

O 

O 


CM 

o 

o 


CM 

,2 

o 


CM 

O 

o 


CM 

O 

O 


0.3 1 5 

Energy (keV) 

Fig. 9.— : The residuals of fits to sample LMXBs with 
the double thermal model with the temperatures of the 
thermal components fixed at values typical of Z sources 
(/cT MC d = 1.7 keV and kT BB =2.5 keV). 

energy band. Using the brightest observation (suz4) of 
4U 1705-44 in ILin etal .l (|201(J ). which has /cTmcd = 1-66 
keV and a BB fraction of 30% in 2-10 keV, we simulated 
200 spectra with Galactic absorption in the direction of 
NGC 3115 and then fitted them with the double ther¬ 
mal model with the column density fixed at the Galactic 
value. We obtained a median of /cTmcd = 1.48t 0 ;y4 keV 
(the error bar corresponds to the 90% confidence level 
in the sense that the upper and lower error bars each 
include 45% of the fits). We also found that 77% and 
46% of the simulated spectra have the best-fitting BB 
fraction in 2-10 keV larger than 30% (the BB fraction 
of suz4) and 50%, respectively. The above simulation 
results indicate that the use of the limited energy band 
could systematically infer a lower disk temperature and 
more BB contribution in 2-10 keV than expected. 

We investigated whether the sources could be fitted 
well if the temperatures of the double thermal compo¬ 
nents are fixed at values typically seen in NS LMXBs. 
We concentrate on six bright sources (the two bright¬ 
est persistent field LMXBs (S25 and S68), the brightest 
transient field LMXB (S35), the two brightest persistent 
LMXBs in GCs (S200 in a blue GC and S212 in a red 
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GC) and the brightest transient GC LMXB (S 8 ); see Ta¬ 
ble [7]). Because they all have near or super-Eddington 
luminosity (Lx > 10 38 erg s^ 1 , brighter than the bright¬ 
est observation (suz4) of 4U 1705-44 in lLin et al.1 (120101) h 
we expect them to be Z sources if they are NS LMXBs. 
Then they should have /cTmcd around 1.7 keV and /cTbe 
around 2.5 keV based on the double thermal model 
(|Lin et al.1 120091 l2012alf . Therefore we also fitted them 
with the temperatures of the thermal components fixed 
at these values. The fit results are given in Table [71 and 
the residuals are shown in Figure[9] The fits overall seem 
fine, but below around 1 keV, we can see a clear system¬ 
atic soft excess in all sources. 

One possible explanation of such a soft excess is that 
the real disk spectrum is broader than the simple MCD 
description. A detailed disk model should carefully cal¬ 
culate radiative transfer through the vertical structure 
of the disk and account for the relativistic smearing. 
This is implemented in the relativistic disk model bh- 
spec, which indicates that the real disk spectrum should 
have excess emission relative to the MC D descr iption 
by more than 10 % at low energies (IDavis et al.1 120051 
Davis fc Hubenvl 120061 : IDavis et al.l [2006: Ku bota et al l 
201011 . Because NSs tend to have a hotter disk, resulting 
in more of the low-energy part of the disk spectrum to 
be observed in the Chandra bandpass than would be the 
case for BHs, the soft excess of the real disk spectrum 
relative to the MCD description is expected to be more 
obvious for NSs than for BHs. Such soft e xcess was not 
seen bv lLin et al.l (I2007L [2009112010L l2012af) because they 
used spectra above 1 keV and the sources that they stud¬ 
ied had much higher absorption than the sources studied 
here. However, bhspec has two limitations, preventing us 
from applying it to our sources. One is that it is a table 
model for the BH accretion disk with the minimum BH 
mass Mbh of 3M 0 , not suitable for the accretion disk 
around NSs. Besides, the model assumes that the inner 
disk radius is at the ISCO, which is not necessarily the 
case for our luminous sources, whose inner disk could 
reach the local Eddington limit and thus be truncated 
outside the ISCO (|Lin et al.ll2009D . 

Given the above possible soft excess problem, we did 
not use the MCD+BB fits to find more BHCs, other than 
those identified in Section 13.31 based on their very soft 
spectra, a characteristic that is not seen in NS LMXBs. 

5. CONCLUSIONS 

We have studied LMXBs detected in NGC 3115 using 
the Megasecond Chandra XVP observations. Including 
three previous observations, the total exposure time is 1.1 
Ms. Thus NGC 3115 is one of Chandra’s best observed 
galaxies. In total we have 136 candidate LMXBs in the 
field and 49 in GCs detected above 2a, with Lx in the 
range of ~10 36 erg s _1 to 10 39 erg s _1 . We calculated 
the long-term variability for all sources and identified 13 
transient candidates, whose long-term variability factors 
are > 5. Excluding these transients, the sources have 


long-term variability overall decreasing with increase in 
luminosity, at least at Lx > 2 x 10 3 ' erg s _1 . 

We carried out simple fits to our sources using single¬ 
component models (a simple PL or a MCD). We found 
that in the Lx versus Tpl and Lx versus /cTmcd plots a 
majority of our sources fall on a narrow track, showing 
harder spectra at higher luminosity below Lx ~ 7 x 10 3 ' 
erg s _1 but relatively constant spectral shape (Tpl ~1.5 
or kTyiCD ~ 1.5 keV) above this luminosity. Because 
our simulations showed that Galactic NS LMXBs in the 
soft state show similar spectral evolution in the Chandra 
bandpass, we identified the track as the NS LMXB soft 
state track and suggested sources below Lx ~ 7 x 10 37 
erg s _1 as atolls in the soft state and sources above this 
luminosity as Z sources. However, the spectra of our 
sources seem to show systematic soft excess rel ative to 
the double thermal (MCD +BB) modeling by ILin et a.l~l 
(12001 [20091 [Mol I20l2al) of Galactic NS LMXBs. One 
explanation for this is that the real disk spectrum has 
excess soft emission relative to the MCD description. 
This is expected from detailed simulations involving care¬ 
ful calculation of radiative transfer through the vertical 
structure of the disk and a c counting for the relativ istic 
smearing (IDavis et al.1120051 IDavis fc Hubenvl 1200611 . 

Ten sources are significantly softer than others at sim¬ 
ilar luminosities and are strong BHCs in the thermal 
state. Five of them are persistent (one in a blue GC), 
and the other five are transient. 

Some special objects were discovered. S198 is the only 
transient BHC whose outburst was covered relatively well 
by our observations. The source displayed clear spectral 
change during the outburst. In the peak, it is consistent 
with a BHC of Mbh ~ 6 Mq in the thermal state with 
fcTiviCD = 0.78 ± 0.07 keV and Lx at 0.12 Eddington 
luminosity. The spectrum during the rise and decay was 
much softer, with kTyiCD = 0.17 ±0.05 keV and the disk 
flux fraction > 63% (the 90% upper limit), and Lx was 
at ~0.008 Eddington luminosity, indicating the possible 
presence of a strong thermal disk at a very low Eddington 
ratio in a BHB. The source could be in the transitional or 
hard state. S109 is a persistent supersoft X-ray source, 
with /cTbb = 86 j ^2 eV and Lbb^oI = O. 8 + 0 I 2 x 10 38 
erg s , and can be explained as due to steady nuclear 
burning on the surface of a WD. S223 is a persistent 
luminous source in a GC with very hard (r = 0.6±0.1) 
X-ray spectra source. It is a candidate NS LMXB with 
a high magnetic field or a high inclination. 
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TABLE 2 

The GC LMXBs and Candidates 


Source 

(i) 

GC 

(2) 

<5xo 

(arcsec) 

(3) 

Sacs 

(mag) 

(4) 

"ACS 

(mag) 

(5) 

<?SCam 

(mag) 

(6) 

^SCam 

(mag) 

(7) 

^SCam 

(mag) 

(8) 

Rh 

(pc) 

(9) 

Vel 
(km s — 
(10) 

4 

A4 

0.06 

20.087 

± 

0.002 

18.786 

± 

0.002 

20.021 

± 

0.003 

19.268 

± 

0.004 

18.968 

± 

0.003 

2.49 

1123 

8 

A238 

0.03 

23.538 

± 

0.021 

21.930 

± 

0.013 

23.766 

± 

0.049 

22.851 

± 

0.045 

22.414 

± 

0.046 

0.94 


10 

A41 

0.19 

21.296 

± 

0.007 

20.250 

± 

0.010 

21.085 

± 

0.107 

20.413 

± 

0.122 

20.274 

± 

0.116 

1.30 

456 

11 

A60 

0.03 

21.896 

± 

0.013 

20.520 

± 

0.015 

21.446 

± 

0.165 

20.704 

± 

0.175 

20.631 

± 

0.183 

1.36 


22 

A171 

0.09 

22.972 

± 

0.018 

21.486 

± 

0.015 

23.067 

± 

0.032 

22.253 

± 

0.033 

21.821 

± 

0.031 

0.77 


23 

A15 

0.03 

20.782 

± 

0.003 

19.473 

± 

0.003 

20.672 

± 

0.005 

19.937 

± 

0.005 

19.644 

± 

0.005 

1.91 

696 

24 

A62 

0.02 

22.071 

± 

0.041 

20.575 

± 

0.042 












70 

A54 

0.02 

21.543 

± 

0.005 

20.429 

± 

0.005 

21.406 

± 

0.008 

20.839 

± 

0.009 

20.633 

± 

0.010 

2.23 

407 

76 

A46 

0.02 

21.858 

± 

0.029 

20.293 

± 

0.028 

20.628 

± 

0.206 

19.748 

± 

0.196 

19.854 

± 

0.239 



96 

All 

0.06 

20.222 

± 

0.004 

19.257 

± 

0.006 










2.59 


101 

A8 

0.05 

20.167 

± 

0.002 

19.013 

± 

0.002 

20.209 

± 

0.004 

19.486 

± 

0.004 

19.217 

± 

0.004 

2.47 


103 

A35 

0.03 

21.619 

± 

0.006 

20.116 

± 

0.005 

21.546 

± 

0.012 

20.699 

± 

0.010 

20.336 

± 

0.011 

2.42 

821 

106 

A7 

0.04 

20.281 

± 

0.002 

18.903 

± 

0.002 

20.263 

± 

0.004 

19.506 

± 

0.004 

19.174 

± 

0.004 

1.76 

697 

113 

A16 

0.12 

20.730 

± 

0.004 

19.503 

± 

0.004 










1.22 


114 

A45 

0.10 

21.612 

± 

0.006 

20.275 

± 

0.005 

21.562 

± 

0.008 

20.824 

± 

0.009 

20.484 

± 

0.009 

2.17 


121 

A10 

0.07 

20.310 

± 

0.003 

19.148 

± 

0.004 

20.138 

± 

0.028 

19.457 

± 

0.030 

19.205 

± 

0.023 

2.42 


129 

A2 

0.13 

19.928 

± 

0.002 

18.752 

± 

0.002 

19.989 

± 

0.003 

19.308 

± 

0.004 

18.992 

± 

0.003 

1.69 


135 

A29 

0.21 

20.981 

± 

0.005 

20.055 

± 

0.006 

21.044 

± 

0.008 

20.412 

± 

0.009 

20.225 

± 

0.009 

1.84 


145 

A57 

0.10 

21.923 

± 

0.006 

20.488 

± 

0.004 

22.024 

± 

0.010 

21.233 

± 

0.010 

20.875 

± 

0.009 

1.08 

798 

150 

A33 

0.13 

21.067 

± 

0.004 

20.084 

± 

0.004 

21.074 

± 

0.007 

20.454 

± 

0.008 

20.252 

± 

0.008 

2.30 

949 

153 

A5 

0.06 

20.033 

± 

0.002 

18.802 

± 

0.002 

20.091 

± 

0.004 

19.360 

± 

0.004 

19.062 

± 

0.003 

1.74 

682 

171 

A14 

0.03 

20.595 

± 

0.003 

19.468 

± 

0.002 

20.623 

± 

0.004 

19.924 

± 

0.005 

19.660 

± 

0.004 

1.57 


183 

A63 

0.06 

21.549 

± 

0.005 

20.585 

± 

0.005 

21.584 

± 

0.008 

20.958 

± 

0.009 

20.745 

± 

0.009 

1.84 

441 

187 

A99 

0.11 

22.531 

± 

0.010 

20.975 

± 

0.006 

22.525 

± 

0.013 

21.679 

± 

0.014 

21.292 

± 

0.013 

2.12 

609 

188 

A13 

0.03 

20.434 

± 

0.003 

19.424 

± 

0.003 

20.394 

± 

0.004 

19.747 

± 

0.005 

19.513 

± 

0.004 

1.24 

426 

192 

A262 

0.03 

23.616 

± 

0.022 

22.202 

± 

0.019 

23.671 

± 

0.034 

22.909 

± 

0.040 

22.530 

± 

0.041 

1.24 


200 

A17 

0.07 

20.862 

± 

0.003 

19.552 

± 

0.003 

20.767 

± 

0.005 

20.003 

± 

0.005 

19.692 

± 

0.005 

2.54 

439 

212 

A102 

0.01 

21.873 

± 

0.006 

21.002 

± 

0.006 

21.838 

± 

0.008 

21.285 

± 

0.010 

21.111 

± 

0.010 

2.56 


213 

A232 

0.15 

22.836 

± 

0.011 

21.887 

± 

0.011 

22.913 

± 

0.018 

22.303 

± 

0.022 

22.103 

± 

0.023 

1.44 

885 

214 

A167 

0.12 

22.901 

± 

0.011 

21.469 

± 

0.008 

22.969 

± 

0.016 

22.178 

± 

0.018 

21.843 

± 

0.017 

1.10 

750 

219 

A339 

0.09 

24.323 

± 

0.036 

23.173 

± 

0.035 

24.557 

± 

0.059 

23.568 

± 

0.058 

23.520 

± 

0.071 



223 

A32 

0.08 

21.179 

± 

0.004 

20.080 

± 

0.004 

21.235 

± 

0.006 

20.548 

± 

0.007 

20.297 

± 

0.007 

1.44 

706 

226 

A24 

0.09 

20.895 

± 

0.003 

19.952 

± 

0.003 

20.913 

± 

0.005 

20.307 

± 

0.006 

20.113 

± 

0.006 

2.37 

809 

229 

A297 

0.04 

23.896 

± 

0.026 

22.597 

± 

0.023 

23.895 

± 

0.035 

23.458 

± 

0.052 

22.999 

± 

0.046 



285 

A249 

0.37 

22.933 

± 

0.013 

22.051 

± 

0.016 

22.935 

± 

0.019 

22.393 

± 

0.024 

22.185 

± 

0.025 

2.54 


299 

A18 

0.15 

20.981 

± 

0.004 

19.641 

± 

0.003 

20.877 

± 

0.005 

20.109 

± 

0.005 

19.805 

± 

0.005 

2.37 

881 

358 

A192 

0.97 

22.516 

± 

0.008 

21.628 

± 

0.008 

22.528 

± 

0.012 

22.044 

± 

0.016 

21.797 

± 

0.015 

4.78 

618 


199 

322 

325 

332 

356 

378 

451 

S624 

S570 

S454 

S578 

S364 

S664 

S638 

0.12 

0.05 

0.64 

0.53 

0.64 

0.59 

0.87 



22.176 ± 0.010 
23.990 ± 0.031 
22.852 ± 0.015 
22.729 ± 0.014 
22.867 ± 0.015 
22.220 ± 0.010 
19.568 ± 0.003 

21.403 ± 0.010 
23.158 ± 0.032 
22.257 ± 0.018 
21.953 ± 0.015 
22.105 ± 0.016 
21.631 ± 0.012 
18.802 ± 0.003 

21.051 ± 0.009 
22.797 ± 0.029 
21.991 ± 0.016 
21.597 ± 0.013 
21.814 ± 0.015 
21.379 ± 0.011 
18.498 ± 0.003 


782 

12 


0.06 

23.377 ± 0.045 

21.832 ± 0.042 




0.25 


53 


0.02 

21.751 ± 0.062 

20.378 ± 0.069 




1.43 


79 


0.03 

21.053 ± 0.033 

19.678 ± 0.033 




1.01 


92 


0.07 

21.693 ± 0.006 

20.306 ± 0.005 

21.738 ± 0.011 

20.959 ± 0.011 

20.608 ± 0.010 

1.66 

238 

Note. 

— Columns: (1) master source unique 

index; (2) GC in 

1 Jennings et al.l ( 20141); (3) the offset between the GC center and 

our X-ray source; 


(4) HST /ACS g-band magnitude; (5) HST /ACS 2 -band magnitude; (6) SCam g- band magni t ude; ( 7) SCam r-band magnitude; (8) SCam i-band 
magnitude; (9) half-light radius; (10) heliocentric velocity, if available, from the IPota et al.1 <120131) catalog. The top group includes 37 LMXBs 
coincident with HST /ACS GCs. The middle group includes 7 LMXBs coincident with GCs detected/covered only in the Subaru/Suprime-Cam 
images, but not in the HST /ACS images. The bottom group includes 4 GC LMXB candidates whose optical matches were not classified as GCs 
by I Jennings et al.l d2014f) but were assumed to be GCs by us (there is another similar source, S65, which is not included in the table because the 
photo metry is not ava il able d ue to being too close to the bright galaxy center; see text for details). Our GC LMXB list is slightly different from 
that in I Jennings et al.l d2014j) , mainly because of our exclusion of very faint X-ray sources (<2cr) and update of calibration in our X-ray analysis. 
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TABLE 3 

The master source catalog 


Source 

CXOU Name 

PU 

Qf/i?25 

-^X,max 

S/N 

^var 

G Lmax 

Type 

(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

(6) 

( 7 ) 

(8) 

( 9 ) 

78 

J100515.4-074254 

0.05 

0.168 

7.41e+37 

20.3 

9.2 

1 

F,BH 

96 

J100516.2-074235 

0.05 

0.230 

4.33e+37 

19.5 

1.7 

2 

GC,BH 

97 

J100514.2-074233 

0.05 

0.279 

5.41e+37 

19.0 

2.6 

7 

F,BH 

100 

J100517.1-074217 

0.07 

0.323 

2.67e+37 

13.3 

2.1 

2 

F,BH 

104 

J100516.5-074207 

0.06 

0.344 

6.40e+37 

19.1 

4.0 

2 

F,BH 

108 

J100518.5-074138 

0.06 

0.518 

1.27e+38 

19.3 

1.3 

0 

F,BH 

179 

J100510.9-074533 

0.18 

1.040 

4.83e+37 

6.2 

19.8 

1 

F,BH,T 

181 

J100510.0-074529 

0.15 

0.936 

1.71e+38 

7.9 

126.9 

3 

F,BH,T 

193 

J100508.7-074443 

0.28 

0.571 

3.99e+37 

2.9 

30.1 

1 

F,BH,T 

198 

J100506.7-074433 

0.08 

0.711 

1.17e+38 

20.2 

33.8 

1 

F,BH,T 

8 

J100517.2-074352 

0.05 

0.913 

1.60e+38 

37.2 

54.8 

2 

GG,T 

25 

J100516.7-074317 

0.04 

0.532 

1.88e+38 

42.6 

1.5 

2 

F 

35 

J100515.8-074312 

0.04 

0.349 

2.10e+38 

37.1 

40.4 

2 

F,T 

68 

J100513.7-074300 

0.03 

0.078 

2.28e+38 

46.7 

1.5 

0 

F 

200 

J100506.0-074428 

0.06 

0.808 

5.71e+38 

69.2 

1.4 

0 

GC 

212 

J100527.3-074316 

0.07 

2.234 

4.85e+38 

57.7 

1.4 

2 

GC 


Note. - This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion 
is shown here, using 10 BHCs (the top group) and 6 bright candidate NS LMXBs (the bottom group), for guidance 
regarding its form and content. Columns: (1) master source unique index; (2) IAU name (following the convention of 
CXOU Jhhmmss.s +/-ddmmss ): (3) 1-a statistical positional uncertainty (in units of arcsec) in each coordinate, based 
on Equation (12) of iKim et al.l (120071) : (4) The ratio of the angular offset from the galaxy center to the elliptical radius 
of the Z ?25 isophotal ellipse in the direction from the galaxy center to the source; (5) 0.5-7 keV maximum luminosity 
(in units of erg s -1 , assuming a source distance of 9.7 Mpc); (6) the signal to noise ratio (the 0.5-7 keV net counts 
divided by the error); (7) long-term variability; (8) the maximum Gregory-Loredo short-term variability index among 
different observations; (9) the source type (field LMXB (“F”), GC LMXB (“GC”), transient (“T”), BHC (“BH”), 
star, AGN, galaxy (“G”), and SSS). 
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TABLE 4 

The Source Counts, Flux, and Hardness Ratio in the Merged Observation 


Source 

Obs 

Expo 

c h 

cl 

C b u 

F h 

K 

K 

hr 3 

hr 3 

HR 3 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

(12) 

78 

0 

1131.0 

418.4 

397.8 

439.2 

3.36e-15 

3.20e-15 

3.53e-15 

-0.19 

-0.24 

-0.14 

96 

0 

1131.0 

385.9 

366.1 

405.9 

3.16e-15 

3.00e-15 

3.32e-15 

-0.23 

-0.28 

-0.18 

97 

0 

1131.0 

363.2 

344.1 

382.5 

2.96e-15 

2.81e-15 

3.12e-15 

-0.34 

-0.39 

-0.29 

100 

0 

1131.0 

181.8 

168.1 

195.6 

1.47e-15 

1.36e-15 

1.58e-15 

-0.53 

-0.60 

-0.47 

104 

0 

1131.0 

366.1 

346.9 

385.5 

2.96e-15 

2.81e-15 

3.12e-15 

-0.19 

-0.24 

-0.13 

108 

0 

443.2 

367.1 

348.0 

386.3 

9.98e-15 

9.46e-15 

1.05e-14 

-0.27 

-0.32 

-0.22 

179 

0 

1131.0 

50.1 

42.0 

58.2 

4.40e-16 

3.70e-16 

5.11e-16 

-0.45 

-0.59 

-0.33 

179 

h 

116.3 

45.2 

38.6 

52.3 

3.51e-15 

3.00e-15 

4.07e-15 

-0.43 

-0.58 

-0.32 

181 

0 

1131.0 

73.3 

64.0 

82.7 

6.79e-16 

5.93e-16 

7.66e-16 

-0.18 

-0.30 

-0.07 

181 

h 

35.8 

74.9 

66.3 

83.6 

1.51e-14 

1.34e-14 

1.69e-14 

-0.06 

-0.18 

0.05 

193 

0 

1061.3 

13.8 

9.0 

19.2 

1.30e-16 

8.50e-17 

1.80e-16 

-0.05 

-0.33 

0.20 

193 

h 

35.8 

17.5 

13.5 

22.1 

3.54e-15 

2.74e-15 

4.46e-15 

-0.14 

-0.41 

0.07 

198 

0 

1131.0 

413.1 

392.6 

433.8 

3.82e-15 

3.63e-15 

4.Ole-15 

-0.19 

-0.25 

-0.14 

198 

h 

414.9 

393.2 

373.5 

413.2 

9.68e-15 

9.19e-15 

1.02e-14 

-0.16 

-0.22 

-0.11 

198 

1 

564.0 

22.5 

17.4 

28.2 

3.98e-16 

3.08e-16 

5.00e-16 

-0.77 

-0.94 

-0.68 

8 

0 

1131.0 

1370.0 

1333.2 

1407.2 

1.10e-14 

1.07e-14 

1.13e-14 

-0.05 

-0.09 

-0.02 

8 

h 

978.9 

1370.8 

1334.0 

1407.9 

1.30e-14 

1.26e-14 

1.33e-14 

-0.07 

-0.10 

-0.03 

25 

0 

1131.0 

1789.1 

1747.1 

1831.5 

1.44e-14 

1.40e-14 

1.47e-14 

-0.09 

-0.12 

-0.06 

35 

0 

1131.0 

1356.0 

1319.5 

1393.0 

1.10e-14 

1.07e-14 

1.13e-14 

-0.10 

-0.13 

-0.06 

35 

h 

978.9 

1144.2 

1110.6 

1178.2 

1.09e-14 

1.06e-14 

1.12e-14 

-0.10 

-0.13 

-0.06 

68 

0 

1131.0 

2186.4 

2139.7 

2233.7 

1.76e-14 

1.72e-14 

1.80e-14 

-0.13 

-0.16 

-0.11 

200 

0 

1131.0 

4716.3 

4648.2 

4785.1 

4.34e-14 

4.28e-14 

4.40e-14 

0.00 

-0.02 

0.02 

212 

0 

1130.8 

3280.4 

3223.5 

3337.8 

3.45e-14 

3.39e-14 

3.51e-14 

-0.04 

-0.06 

-0.01 


Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion 
is shown here, using 10 BHCs (the top group) and 6 bright candidate NS LMXBs (the bottom group), for guidance 
regarding its form and content. Columns: (1) master source unique index; (2) observation (“0” refers the combination 
of all available observations; “h” refers to the combination of the high-state observations for the 13 transients, as 
shown in Figure [3] a-d; and “1” refers to the combination of the low-state observations for source 198, as shown in 
Figure l3cl): (3) exposure time, in units of ks; (4-6) the net counts in the broad band (0.5-7.0 keV) and the lower and 
upper limits; (7-9) the energy flux in the broad band and the lower and upper limits, in units of erg s^ 1 cm -2 ; and 
(10-12) the hardness ratio using the fluxes in the 0.5-2.0 keV and 2.0-7.0 keV energy bands and the lower and upper 
limits. All limits are at the 68% confidence level. 
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TABLE 5 

The Source Counts and Flux in Individual Observations 



ObsID 

c b 

Cl 


F h 


F b 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

78 

13820 

86.9 

77.6 

96.4 

4.37e—15 

3.90e—15 

4.84e—15 

96 

13820 

49.0 

42.2 

56.5 

2.47e—15 

2.13e—15 

2.85e—15 

97 

13820 

56.0 

48.5 

63.5 

2.82e—15 

2.44e—15 

3.20e—15 

100 

13820 

32.1 

26.6 

38.1 

1.60e—15 

1.33e—15 

1.91e—15 

104 

13820 

38.6 

32.5 

45.2 

1.93e—15 

1.63e—15 

2.26e—15 

179 

13820 

0.7 

0.0 

2.8 

3.78e—17 

0 

1.46e—16 

181 

13820 

0.0 

0.0 

1.1 

0 

0 

5.97e—17 

193 

13820 

0.0 

0.0 

1.1 

0 

0 

5.87e—17 

198 

13820 

174.9 

161.7 

188.2 

1.04e—14 

9.60e—15 

1.12e—14 

8 

13820 

237.2 

221.9 

252.7 

1.19e—14 

1.12e—14 

1.27e—14 

25 

13820 

294.4 

277.4 

311.7 

1.48e—14 

1.39e—14 

1.56e—14 

35 

13820 

200.9 

186.8 

215.1 

1.09e—14 

1.Ole—14 

1.17e—14 

68 

13820 

366.4 

347.2 

385.7 

1.84e—14 

1.75e—14 

1.94e—14 

200 

13820 

707.1 

680.7 

733.7 

4.23e—14 

4.08e—14 

4.39e—14 

212 


461.3 

439.9 

482.8 

3.17e—14 

3.02e—14 

3.31e—14 


Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion 
is shown here, using 10 BHCs (the top group) and 6 bright candidate NS LMXBs (the bottom group) in observation 
13820, for guidance regarding its form and content. Columns: (1) master source unique index; (2) observation ID; 
(3-5) the net counts in the broad band (0.5-7.0 keV) and the lower and upper limits; and (6-8) the energy flux in the 
broad band and the lower and upper limits, in units of erg s -1 cm -2 . All limits are at the 68% confidence level. 
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TABLE 6 

Spectral fit results 





the PL Model 



the MCD Model 


Source 

Obs 

TVS 

r 

JNorm 

L 

7V+ 

kT 


L 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

78 

0 

18.0+1 

2.84+; 3 2 

1 *31 +0-39 
± * ol -0.28 

q 1 +0.5 
°- 1 -0.3 

0.0+ 21 

0.641+ 7 

.39.09++ 

3-0+1 

96 

0 

14.4+1 

2 7Q +0 - 35 
z - ‘^-0.33 

1 1 O+0-32 
1 - lz -0.24 

3-0+1 

0.0+ 14 

0 62 +0 ' 07 
u - uz —0.06 

40.651+7 

9 o+0.3 
^•°-0.4 

97 

0 

16-7 ±?;g 

3-36++ 

1 nr+0.42 
i.zO-o.30 

2 - 8+1 

0.0+ 14 

0-45+1 3 

76.14++8 

2 6 +0 ' 3 
z -°-0.2 

100 

0 

o.o+ 3 - 5 

3.08+1° 

0.38+;°® 

9 q+0.3 
Z ’°-0.4 

o.o+ 10 

0-291+4 

152.16+°+ 

1 3+°- 3 
± *°-0.2 

104 

0 

16-5+1 

2 7fi +0 ' 37 

z *' °-0.34 

i-09±g;i| 

2 - 8+1 

o.o+ 2 - 3 

0.641+® 

36.70+; 34 

2 6 +0,2 
z - u -0.2 

108 

0 

24.9+1 

O 24+0.40 
O.Z^i— o 37 

4.95++ 

8-7+1 

0 0+3.9 

0.591+7 

79.87I43I0 

8.7+1 

179 

h 

qi -i +26.5 
01.1—26.1 

4 21 +0-79 

+.Z1-1 44 

2-26tl.ll 

9 Q + l-3 
^•°- 0.8 

0.0+ 201 

n qq+ 013 

u.ob;_o.i4 

1 1 1 32+ 279 - 70 
iii.oz_ 51 07 

3.0+1 

181 

h 

10 l+ 13 - 0 

2 96 +0 ' 84 
z. JU-0.69 

4-58++ 

14.0+1 

0.0+ 3 ' 5 

n 43+0.11 

u.+o_o.08 

171 ni+ 92 - 25 

1 1 i.ui_ 59 61 

12 . 4+1 

193 

h 

9-0 ++ 

2 46 +2-54 
z * w -1.03 

0-94++ 

q q+2.7 
°*°- 1.8 

o.o + 310 

0.571+ 3 

4q 1 K+201.12 
^y.io_ 31 44 

2.8+1 

198 

h 

20.3+3 

2 3q+°- 34 
z -°^-0.32 

3 63+ 1 - 12 
O.DO-o.84 

q Q + 1 - 4 
y.O_i 2 

0.0+ 3 ' 5 

0-80+1° 

43.19++ 6 

9.0+1 

198 

1 

20.3(f) 

5.00 o.8o 

n i q+0-08 

0.Io 0.06 

0-5+1 

0.0(f) 

n 1 7+0.07 
u - 4 ‘ -0.04 

345.30I235.01 

0 6 +0,4 
U.O-o 2 

8 

h 


i ry+u.u 
± - 0 ' -0.13 

9 £7+U-3t> 

-0-32 

18-0+1 

0.0 +ut> 

1 47 +u.ics 
-Q.ii 

17.12++ 

16.21+ 

25 

0 

q q+3.5 

°* z -3.2 

1 43+°- 11 
± ’™— 0.11 

2 34+ 0 - 29 
z -°^-0.27 

21.31+ 

o.o + 0 - 6 

1 62 +0 ' 14 
i ,uz —0.10 

15.31+;°° 

18.9+1 

35 

h 

3-6±|;i 

1 34+0-14 
1 -°^-0.14 

2 03+°- 31 
z .UO_ 0 .27 

15-1+1 

0 o+°- 6 

1 43+0.14 
4-^0—0.12 

15.91+; 7 ° 

13.4+ 3 

68 

0 

i-s+1 

1 48 +0 - 10 
± -^°-0.09 

3 in+0-33 
°- lu -0.27 

25.21+ 

o.o + 0 - 4 

1 31 +0-H 
-L.01_o.io 

19.08++ 

99 9+2-1 
ZZ ’ Z — 1.0 

200 

0 

5-7+1 

1 4fi+°- 07 
i ^°-0.07 

o i ,4+0.59 
0.14-0.55 

63.71+ 

o.o +0 - 4 

1 f)1+°- 08 
1 - Oi -0.02 

27.02+;“ 

57.0+1 

212 

0 

q 0+3.7 
y '°-3.6 

1 62 +0 ' 09 
1 - uz -0.09 

7 0Q +0 - 74 
' - uy -0.61 

45.41+ 

0 0+o.e 

1 cq+0.09 
i.OO_o.o 2 

25-541°+ 

42.1+1 


Note. - This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion 
is shown here, using 10 BHCs (the top group) and 6 bright candidate NS LMXBs (the bottom group), for guidance 
regarding its form and content. Columns: (1) master source unique index; (2) the observation (“0” refers the combi¬ 
nation of all available observations; “h” refers to the combination of the high-state observations for the 13 transients, 
as shown in Figure El)a)-(d); and “1” refers to the combination of the low-state observations for source 198, as shown 
in Figure [dfc)); (3)-(6) the intrinsic column density (in units of 10 20 cm -2 , constrained to be < 10 23 cm -2 ), pho¬ 
ton index (constrained to be < 5.0), normalization (in units of 10 -6 photons keV -1 cm~ 2 s _1 at 1 keV), and 0.3-8 
keV luminosity (in units of 10 37 erg s _1 , assuming a source distance of 9.7 Mpc) corrected for Galactic absorption 
from the fit with the PL model; (7)-(10) the intrinsic column density, the maximum disk temperature (in units of 
keV, constrained to be < 4 keV), the apparent inner disk radius 7?mcd (in units of km) from the normalization 
Amcd = ((i?MCD/km)/(U/10kpc)) 2 cos0, where D is the source distance 9.7 Mpc and 9 is the disk inclination as¬ 
sumed to be 60°, and 0.3-8 keV luminosity corrected for Galactic absorption from the fit with the MCD model. All 
fits used the C statistic. All errors are at the 90% confidence level. For S198, the fit to the low-state spectrum has the 
column density fixed (marked with “f”) at the value obtained from the fit to the high-state spectrum. 
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TABLE 7 

Spectral fit results of some special sources 


Source 

model 

n H 




Other Parameters 



xl/v 





( 10 20 cm -2 ) 








1 0 36 

erg s 1 

(i) 

( 2 ) 

(3) 




(4) 



(5) 

( 6 ) 

(7) 

100 

MCD+PL 

0 o + 7 - 2 

fcT MCD = 014 io.04 

keV 

%CD 

= 66 O+I 3® 3 km, r PL = 2.5(/), AT PL 

= (0.26 

CO 

1 

0 

X 

N N 

OO 

do 

+ 1 


22+ 4 

22 + 1 5 

109 

BB 

0 . 7 + 10.6 

feT BB = 0.086+“ ""' keV, fi BB = 3376+“ 2 ^ km 


37+0 

4l" t Yy 

8 

MCD+BB 

0 . 0+ 3 ' z 

feT MCD - ° - 67 —0.20 

keV 

%CD 

= 56^2° km, fcTgg = 1.43^J '° keV 

h bb 

14 - 6 -ll '3 km 

0.93(55) 

186+lg 

1 88 —12 

25 

MCD+BB 

0 . 0+ 11 ' 2 

fcT MCD = °- 46 ioi 8 

keV, 

%CD 

= 100 +|? 8 km, fcT BB = 1.11 + 3 ' 32 keV, H RB 

= 26 -7±io.s k ” 

1.17(72) 

1841“ 

184+g 

35 

MCD+BB 

0 . 0+ 5 - 8 

fc %CD = 

keV 

r MCD 

= 77+ 88 km, feT BB = 1 . 22 +°^ keV 

fi BB 

18.4t 3 Q 4 0 km 

0.74(45) 

“Til? 

137t! 2 n 

68 

MCD+BB 

1 7+5.6 

‘%CD = 

keV, -RmCD 

= 154+1 28 km, fcT BB = 1.24+ 8 ' 2 ! keV, H BB 

= 24 , 1 ±5.9 km 

0.82(84) 

23 +J 4 

2 39t 22 

200 

MCD+BB 

0 . 0+ 1 - 0 

fcT MCD = 0-62 i o. 15 

keV, 

fi MCD 

= i06±™ km . fcT BB = i.ie+g-?J k «v, f? BB 

= 4 1' 5 ±2T3 km 

0.96(165) 

556±l\ 

336+21 

212 

MCD+BB 

o.o + 5 - 8 

feT MCD = °- 56 iS'.l5 

keV, %CD 

= 113± 4 ° 0 km, fcT BB = 1.15+q' 2 3 keV, il BB 

= 36.0+J5'4 km 

0.91(123) 

416+ 16 
4ie> —20 

416+30 
41t) —16 

8 

MCD+BB 

0 , 0+ U ' 6 

feT MCD = 1-70(f) 

keV 

%CD 

= 13.1^J| km, fcT BB = 2.50(f) keV, 

R BB = 

0 . 0 + ie km 

1.37(57) 

165+^. 

165+y 

25 

MCD+BB 

0 . 0+ 06 

WmCD = 1-70(f) 

keV 

%CD 

= 13.8 + 8, 4 km, fcT BB = 2.50(f) keV, 

fi BB = 

0.1 + 2 ' 6 km 

1.38(75) 

183± 41 

183+y 1 

35 

MCD+BB 

0 . 0+ 06 

*7 MC D = 1 ■ 70 (f) 

keV 

h mcd 

= 12.0+q | km, fcT BB = 2.50(f) keV, 

fi BB = 

0 . 0+ 1 ' 5 km 

1.32(47) 

137+ 8 

137+3 

68 

MCD+BB 

0 . 0+ 04 

feT MCD = 170(f) 

keV 

%CD 

= 15.2+°'J km, fcT BB = 2.50(f) keV, 

fi BB = 

0.1 + 2 ' 2 km 

1.73(86) 

220 +| 

220 +| 

200 

MCD+BB 

0 . 0+ 04 

WmcD = 1.70(f) 

keV 

fi MCD 

= 24.3 + Q 4 km, fcT BB = 2.50(f) keV, 

r BB = 

0.1 + 2 ' 4 km 

1.05(167) 

565^^3 

365±i| 

212 

MCD+BB 

0 . 0+ 04 

17 MC D = 170(f) 

keV 

%CD 

= 21.0+J5 ' 4 km ’ fcT BB = 2.50(f) keV, 

fi BB = 

0 . 1+ 1 ' 8 km 

1.14(125) 

422+i 2 

422 +; 2 


Note. — Columns: (1) master source unique index; (2) the spectral model; (3) intrinsic column density (4) other spectral parameters (Rgg is the apparent source radius 
from the BB normalization iV BB = ((Hgg/km )/(D / lOkpc)) 2 , and refer to Table [6l for the meanings of other parameters); (5) the reduced x^ and the degrees of freedom, 
(only for fits to spectra binned to have a minimum of 20 counts and using the x^ statistic), (7)—(8) the absorbed (but corrected for Galactic absorption) and unabsorbed 
0.3—8 keV luminosity (assuming isotropic emission or a disk inclination of 60°), respectively. All errors are at the 90% confidence level. 














